Intrinsic dissipation in nanomechanical resonators due to phonon tunneling 
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State of the art nanomechanical resonators present quality factors Q ~ 10"^ — 10^, which are 
much lower than those that can be naively extrapolated from the behavior of micromechanical 
resonators. We analyze the dissipation mechanism that arises in nanomechanical beam-structures 
due to the tunneling of mesoscopic phonons between the beam and its supports (known as clamping 
losses). We derive the environmental force spectral density that determines the quantum Brownian 
motion of a given resonance. Our treatment is valid for low frequencies and provides the leading 
contribution in the aspect ratio. This yields fundamental limits for the Q-values which are described 
by simple scaling laws and are relevant for state of the art experimental structures. In this context, 
for resonant frequencies in the 0.1 — IGHz range, while this dissipation mechanism can limit flexural 
resonators it is found to be negligible for torsional ones. In the case of structureless 3D supports the 
corresponding environmental spectral densities are Ohmic for flexural resonators and super-Ohmic 
for torsional ones, while for 2D slab supports they yield 1// noise. Furthermore analogous results 
are established for the case of suspended semiconducting single-walled carbon nanotubes. Finally, 
we provide a general expression for the spectral density that allows to extend our treatment to 
other geometries and illustrate its use by applying it to a microtoroid. Our analysis is relevant for 
applications in high precision measurements and for the prospects of probing quantum effects in a 
macroscopic mechanical degree of freedom. 
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I. INTRODUCTION 

Nanomec hanic al resonators offer a host of novel 
application^iJ^El \^ high precision measurements^ and 
may provide a new arena for probing fundamental as- 
pects of quantum physic^. A prominent example of such 
novel applications is scaling magnetic resonance force 
microscopy^ (MRFM) down to the level of single-spin de- 
tection. This degree of sensitivity would allow, for exam- 
ple, three-dimensional imaging of individual biomolecules 
with atomic-scale resolution. In these microscopes a 
magnetic particle!^ mounted on a cantilever interacts with 
the nuclear or electron spins in the sample via the mag- 
netic dipole force. Hence as in many other applications 
the relevant mechanical resonator has the beam geometry 
on which we will mainly focus in this paper™. Beyond mi- 
croscopy, an interesting issue is whether in the nanoscale 
regime such a device could be used to probe the quantum 
state of a single spin^°. 

A single spin 1/2 coupled to a harmonic oscillator, 
namely the Jaynes-Cummings model, constitutes a fun- 
damental system in quantum optics^^. Physical realiza- 
tions of this system have proved invaluable in probing and 
understanding quantum phenomena^^ These investi- 
gations have ranged from studying the quantum-classical 
interface to proof-of-principle demonstrations of the ba- 
sic building blocks of quantum information processing. 
These developments were enabled by the advent of high- 
er resonators for optical and microwave photons and of 
conservative harmonic traps for atoms. In these systems 
the relevant harmonic oscillator is furnished, respectively, 
by a single normal mode of the electromagnetic field or 



the atomic motiorP^. On the other hand state of the 
art semiconductor nanostructures can support mechan- 
ical resonances with quality factors Q ~ 10^ — 10^ and 
frequencies approaching the GHz regimd^. For some of 
these structures measurements of mechanical displace- 
ments with a sensitivity approaching the quantum limit 
have been achieved by expl oiting capacitive coupling to 
a single electron transistoil ^^l^^l These developments 
suggest the possibility of realizing a quantum vhonon- 
ics realm in complete analogy to quantum optic^^EH. In 
this case the relevant harmonic oscillator would be fur- 
nished instead by a mechanical resonance and the role 
of the pscudospin co uld b e played by a capacitively cou- 
pled Cooper pair boxE^E^l or an excitonic transition of an 
embedded self-assembled quantum dot!^. Along these 
lines one could envisage observing quantum jumps due 
to the discrete nature of phonons and realizing quantum 
state engineering of non-classical states of motion^^. In 
addition to semiconductor planar heterostructure real- 
izations there are other promising possibilities like sus- 
pended single walled carbon nanotubes^^ ^'^ (SWNT), 
nanowires^Sl, and single crystal diamond beams with em- 
bedded nitrogen vacancy color center^^l Another alter- 
native for furnishing the non-linearity needed to induce 
non-classical behavior in mesoscopic mechanical oscilla- 
tors are optomechanical schemes in which the resonator 
couples via radiation pressure to an optical cavitj^^^l^ll 
This venue has recently witnessed significant experimen- 
tal progres^2ll towards achieving ground state cooling 
which is highly desirable to enable quantum effecta^Sl 
Finally, for resonators with sufficiently small effective 
masses and high bending rigidity, yet another alternative 
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that has been considered is to use the Euler instabilitj!^. 

The reahzation of all of the aforementioned applica- 
tions of nanomechanical systems hinges on understand- 
ing and controlling the intrinsic dissipation and noise 
mechanisms'^^ that limit their coherent dynamics. For 
centimeter-scale semiconductor micromechanical systems 
Q ~ 10® have been measured at low temperature^. 
On the other hand when these devices are shrunk to 
the nanometer-scale these values decrease dramatically 
to Q '--^ 10^ — 10^. Early work suggested that the in- 
crease in surface-to-volume ratio combined with surface 
effects might be invoked as a plausible explanation of this 
phenomenon^. Subsequently it was realized that elastic 
wave radiation into the supports — the so-called clamp- 
ing losses^ — could play an important role leading to 
non-trivial scaling laws for the Q-values with the aspect 
ratio that are intimately related to the low frequency 
behavior of the corresponding transmission coefficients 
(cf. ~ 



Sec. 



IV]and Refs. 33,34 35,^). 
If we consider the environment of the mechanical res- 
onator responsible for its dissipation, the Q- value is de- 
termined by the pole of the ensuing modified propagator 
for the resonator's normal coordinate that corresponds 
to its resonan t freq uency ujh. If one adopts a Caldeira- 
Leggett modeW^^ — i.e. the environment is assumed 
to consist of a thermal ensemble of harmonic oscillators 
— with a linear coupling to the environment, this prop- 
agator can be obtained exactly from the environmental 
force spectral density and the Q- value is found to be 
temperature independent. Thus, the quantum Brownian 
motion of the normal coordinate Xji associated with a 
given resonance is characterized by the following gener- 
alized equation of motion 

+ fdt' 7(i - t')XR{t') + ulXuit) = m (1) 

"'0 

and is completely determined by the function I{lo) = 
Lu J^^df 7(i)e"^'/2wii where 7(t) is the symmetric 

dissipation-kerneP' [cf. Eq. ^ and Sec. |ll] and i{t) 



corresponds to the environmental noise. When the stan- 
dard Markov approximation is warranted j(t) ~ S{t) and 
the Q- value is determined by 



1 

Q 
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Within this approximation the value I{liJr) is the only 
relevant information about the spectral density, Eq. ([T]) 
can be interpreted as a quantum Langevin equation, and 
we have the standard relation 



{E) 



(3) 



where E is the total energy stored in the resonator's de- 
gree of freedom. This is normally valid for sufficiently 
high Q and low temperature if I{lo) is smooth enough, 
or for high temperatures if I{lo) scales linearly with 
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FIG. 1: Schematic diagram of suspended mechanical beam 
structures presenting abrupt junctions with the supports. 
Throughout the paper the z-axis corresponds to the beam's 
axis, while the x and y axes are the principal a:xes of inertia 
of the beam's cross section at the left junction (the origin is 
set at the center of mass of the latter). 



frequencjEHSSIini (classical Langevin). However in many 
instances it is desirable to go beyond the Markov ap- 
proximation and more precise knowledge about /(w) is 
needed. Examples of these are the cases discussed above 
in which the resonator interacts with a single spin 1 /2 
system if the coupling is relatively strong and cases in 
which the Q is not high enough or the temperature low 
enough. Currently there exists no microscopic derivation 
of Eq. Q and of the underlying Caldeira-Leggett Hamil- 
tonian and corresponding spectral density , for which 
customarily an ohmic dependence is assumed. 

For nanometer-sized suspended monocrystalline beam 
structures at low temperatures, the relevant phonon 
mean free path can become larger than the beam's 
length. In the case of an insulating system, a promi- 
nent consequence of the ballistic regime that results 
is the quantization of thermal conduct ancd^ in units 
of Gth = TT^k%T/3h, which has been demonstrated 
experimentally "*^-. In this low temperature, low frequency 
regime and for small deflections, anharmonicity becomes 
irrelevant for the analysis of dissipative eSect^^^^^. 
Thus, in this paper we analyze the ideal limit, that en- 
sues for the low-lying resonances of an insulating beam 
close to equilibrium at low enough temperatures, in which 
phonon tunneling between the beam and its supports is the 
only source of thermal noise and dissipation. To this ef- 
fect the vibrations of the whole structure are described 
by a purely harmonic Hamiltonian. Its normal modes, 
discussed in the following section, form a continuum. As 
will become clear below, this feature is inherited from the 
supports that thus provide thermal phonon reservoirs. 
When the systems under consideration deviate from the 
aforementioned ideal scenario (e.g. the vacuum is imper- 
fect, the beam's material is amorphous or surface effects 
are relevant) there will be other contributions to the me- 
chanical damping that will add incoherently'^. However it 
should be stressed that in all such cases our treatment is 
valid for the contribution to the dissipation arising from 
the vibrational degrees of freedom and thus our results 
will provide an upper bound for the Q-values. 

Beam nanostructures for which high Q-values have 
been measured normally involve abrupt junctions with 
the supports. This can be intuitively understood in terms 
of impedance mismatch for the propagation of elastic 
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waves. We will focus mainly on this type of structures, as 
depicted in Fig. [l] and consider phonon frequencies that 
correspond to fc ~ 27r/L ^ 27r/d, where k is the wavevec- 
tor along z inside the beam, L the beam's length, and 
d the typical dimension of the cross section (cf. Fig. [l] 
and App. |b]). It is clear that for this beam geometry 
the reciprocal of the aspect ratio d/L provides a natu- 
ral small parameter that will underpin our analysis of 
the non-trivial matching at the abrupt junctions. Given 
that we are interested in phonons with wavelengths much 
longer than the lattice constant, the continuum limit 
constitutes a good approximation. This yields for the 
supports 3D isotropic elasticity^ '^^ , if one further ne- 
glects the possible anisotropy — or 2D "thin plate elas- 
ticity" for slab support s of t hickness much smaller than 
the phonon wavelengtlJ22llIl. If we consider the beam 
there are basically two scenarios: (i) all the charac- 
teristic dimensions of the cross section are much larger 
than the lattice constant (e.g. semiconductor planar het- 
erostructures) or (ii) there are characteristic dimensions 
comparable to the lattice constant (e.g. semiconducting 
SWNTs and small-radius nanowires). In case (i) one can 
start from 3D isotropic elasticity and obtain the effective 
theory valid for phonon wavel engths much larger than d 
known as "thin rod elasticity"I^B^ (TRE). The latter 
is completely determined by the linear mass density /i;,, 
the mean axial moment of inertia {r'^)s7 and the exten- 
sional, torsional and bending rigidities; and leads to the 
following dispersion relations for the four low frequency 
branches 



quantum dissipative dynamics of a "macroscopic" me- 
chanical degree of freedom can be g iven, w hich reduces 
exactly to a Caldeira-Leggett modePSMMSJ, xhe latter 
is essentially a derivation of the basic Hamiltonian for 
a lossy ID phonon-cavity and brings together aspects of 
the analogous problem in quantum opticj^ and the bal- 
listic transport in a mesoscopic wire^. In this respect it 
should be noted that the Markov approximation, which is 
customarily used for the Brownian motion of mechanical 
resonators, is not assumed in our derivation but instead 
its range of validity emerges from the behavior of the 
spectral densities I{uj) that we calculate and that allow 
(if necessary) to go beyond the Markovian regime — for 
these aspects the reader is referred to the vast literature 
on quantum dissipation^^. Thus, the validity (for large 
aspect ratios) of Eq. Q , that will be used to determine 
the Q- values, is established rather than postulated. 



We will first present our main results (cf. Subsec. lAl 
and subsequently give their derivation in the following 
Sections. For the sake of clarity we will mainly focus on a 
specific model with maximally symmetric 3D structure- 
less supports and the beam having the aforementioned 
property (z) [cf. Subsec. II A , and outline how the treat 



ment can be extended (cf. App. [C| by referral to this 
concrete realization. We will discuss both the bridge and 
cantilever geometries (cf. Fig. [T). Finally, we will also 
outline how this approach can be generalized to other 
geometries and illustrate this further extension for the 
experimentally relevant case of a microtoroicP^'^. 
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Here /3 is the branch index: c for compression, t for tor- 
sion, V for vertical bending and h for horizontal bend- 
ing [the prefactors C/s are given in App. Ib] Eq. (B9l]. 
Though for propagation inside the beam this effective 
theory will suffice, resort to the underlying 3D "micro- 
scopic" theory will allow for a rigorous treatment of the 
matching at the junction (cf. Sec. IV I — an analogous 
procedure is feasible for a thin plate geometry ■^'^ with the 
"microscopic" theory furnished by 2D thin plate elastic- 
ity. In case (ii) it is well understood that the effective 
long wavelength theory is formally equivalent to TRE 
but the standard "bulk" relations between the rigidities 
are no longer warrantecP^. Whence heuristic consider- 
ations allow to extend our treatment of the non-trivial 
matching at the abrupt junctions to this case with suit- 
ably redefined cp — when considering this extension we 
will focus mainly on SWNTs. Thus, our treatment will 
amount to a derivation of the leading contribution in the 
aspect ratio to the environmental spectral densities I{uj) 
associated with phonon tunneling induced noise starting 
from the underlying microscopic lattice Hamiltonian. 

The relevance of the ensuing results is twofold: {i) they 
furnish a very general understanding of clamping losses 
determining the corresponding Q- values for a wide range 
of experimentally relevant structures, and (ii) they pro- 
vide an instance for which a microscopic derivation of the 



Environmental force spectral densities and 
Q-values for each resonance 



To study the dissipation induced by the coupling to 
the supports it proves useful to introduce the concept of 
an "effective environmental density of states" for each 
branch given by pf3{uj). These functions will be defined 
in Sec. |IV[ There we will find that for 3D supports they 
bear simple relations with properties of the decoupled 
support — i.e. subject to free boundary conditions — 
that are closely related to its density of states (DOS); 
namely, the displacement and angle (twist) vacuum spec- 
trum at the junction^^. The cornerstone of our analysis 
will be furnished by the following novel pair of funda- 
mental relations, each of which completely specifies the 
environmental force spectral density for a given resonance 
to lowest order in the reciprocal of the aspect ratio d/L: 
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General Relations (3D) Monolithic structure {a 



Semiconducting SWNT {a^ = \) 
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TABLE I: First and second columns: general formulas for the environmental force spectral densities /„(a;) and the Q- values 
Qn corresponding to the different resonances of a beam suspended from structureless 3D support(s) — n = 0, 1, . . . labels the 
harmonics for each of the four branches. These general formulas are specialized for two cases: (third column) a monolithic 
structure with rectangular beam cross section of thickness t and width w, and (fifth column) a suspended semiconducting 
SWNT of radius R. The nanotube is modeled as a cylindrical shell of effective thickness h and Poisson raticP^ ab = 0.19. 
All formulas are valid for both the cantilever {5 = 1) and bridge {S — 2) geometries. The dispersion-relation prefactors 
are given in Eq. (B9l, fib is the linear mass density of the beam (for the nanotube p,b = 2tiR(jg, where gq is the surface 
density of graphene), ps and Ct are, respectively, the mass density and transverse speed of sound for the supports' material, 
and (r^) = Iz/ S is the beam's mean axial moment of inertia. In the case of torsion (second row), for the monolithic structure 
(third and fourth columns), we specialize for t ^ w so that ^1 ~ 1 and the torsional rigiditjBS reads C ~ i5t,wt'^/6(l + ah) 
[for a cylindrical shell we have instead C — EtirhR^ / {1 + ab)] — Eb [Eg] is the Youn g's modulus for the material of the beam 



[support(s)]. The dimensionless displacements and angles up, given in Subsec. IVC (cf. Fig. |3|, take the following values for 
^ — 1 /Q r^, — /"I o,^ ^ ^ \ — 1 //il. ^. f^ I A\ — n 1 Q ^. ^^ I a ^, \ — i !^^^ i n m o^,2~^ .^v^^ ^, . /i I a\ — n io 



0.13, Ut(l/4,7, 



) = l/127r + O.OI97J7 and ?i„/h(l/4) = 0.12, where 



1/3 [a EE (1 - 2a,)/2(l - a,) = 1/4]: Mc(1/4) 

7z = [ly — and we have defined ut(a,'yz) = Ut'^' + ul^\a)'yi. The resulting typical values for the Qn (fourth and sixth 

columns) correspond to the lowest lying resonances of bridge geometries (wavevectors" fco,c/t = ^f/i and kg^^/f^ ~ I.SItt/L) 
with a different length for each branch (-L/3) chosen so that the comparison is made for equal frequencies'': 150MHz for the 
monolithic structure (parameters: w = 100 nm and t = 20 nm) and IGHz for the nanotube [parameters: /i = 0.66A(cf. Ref. [52]), 
crG = 7.7xlO~^Kg/m^, £'i, = lTPa (cf. Ref. [23), and i?=lnm; we assume Si supports: = 2.3 x lO^Kg/m^, Es = ll2GI'&]. 
''For the cantilever geometry these would be instead fcp ^./j = ■k/2L and k^^^if^ ^ 0.6Qtt/L. 

''For the monolithic case (nanotube): Lc = 23p (2. 3/^), Lt = 4.7/x (1.5^*), Ly=0.98p. (0.11/i), Lh = 2Afi (O.ll^t). 



Here n = 0, 1, . . . labels the harmonics for each branch 
(i.e. LUn — > uJn,/3), (5 = 1, 2 is the number of supports, 

[1 for/3 = c,t 

and the bare resonant wavevectors fc„_^ = aj^^(aj„_^) are 
determined by the TRE solutions for clamped-clamped 
(clamped-free) boundary conditions in the case of the 
bridge (cantilever) geometrjJ^. The first relation ex- 
presses the spectral density in terms of pp{i^)- The sec- 
ond relation involves the transmission coefficient Tfj^uj) 
from a semi-infinite beam into the support at a single 
junction which can be interpreted as the probability for 
an incident phonon with frequency to to tunnel into the 
support. As will be discussed further in the next Sec- 
tion and derived rigorously in Subsec. |IV A[ for phonon 
frequencies a; — > the beam and the supports become 
decoupled so that the junction plays a role analogous to 



a tunnel barrier. Whence relation (|6| allows to interpret 
the losses to the supports in terms of phonon tunnel- 
ing. More precisely, it provides a rigorous footing for the 
heuristic formula to describe this dissipation mechanism 
set forth by Cross and Lifshitz in Ref. [35] based upon 
Eq. ([3|, namely, 

2L dk^ ""'ujpikn) 2 knL ' ^ ' 

where the implicit dimensionless prefactor is expected 
to be of order unity and have the latter as its limit for 
n — > cx) — henceforth we drop the branch index of the 

km ^n- 

The above approximation emerges from considering a 
phonon wavepacket that bounces back and forth between 
both ends of the beam. Naturally, it should be adequate 
for large n and provides a simple intuitive description of 
this dissipation mechanism in terms of phonon tunneling 
at the junctions. Here we have added a factor of S/2, 
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where 5 is the number of supports, to account for both 
the bridge and cantilever geometries. Equation (|8| can be 
obtained from Eqs. ^ and (|6| which allow in addition to 
determine that the corresponding dimensionless prefactor 
is given by Cn,f3 — i.e. for the non-dispersive branches 
it is exactly unity for all n and only deviates from unity 
for the low-lying bending resonances. 

The general relation between the transmission coeffi- 
cient at a single junction and the effective environmental 
DOS concomitant to the pair of relations ([s]), Q will 
be proved in Sec. IV where expressions for the functions 
T/3(a;) and P0{ lo) are derived explicitly for 3D supports 
[cf. Eqs. ([5§-@]. The latter together with Eqs. ^ 



and ^ allow us to obtain formulas for the Q- values of all 
the low-frequency resonances that are given in the second 
column of Table |l] The first column gives the resulting 
expression for the environmental force spectral density 
which turns out to be ohmic in all cases except for the 
torsional resonances. The third and fifth columns give 
examples of particular experimental relevance. Namely, 
a monolithic structure with rectangular cross section for 
the beam of width w and thickness t and a semiconduct- 
ing suspended SWNT of radius R — for the nanotube 
we use constants that correspon d to t he "continuum" 
shell approximation for the rigiditie^SMH^ r^j^^ validity of 
these results, which are adequate for low frequencies, will 
be borne out in full detail in the following Sections as we 
derive the Eqs. ([S]) and ([6]). While the latter will hold 
in all instances where any characteristic dimension of the 
supports is either much larger than L or at most of order 
ds < d, Table [l] focuses on the case in which the limit 
ds yields a 3D support with no characteristic di- 
mension (e.g. an elastic half-space). Within this context 
the second column is general with an appropriate defi- 
nition of the dimensionless displacements and angles u 
that only depend on the supports' material Poisson ratio 
(cf. Sec. IV and App. [C|. On the other hand the specific 
u used in the examples — which will be calculated in 
Sec. |IV| — correspond to the maximally symmetric case 
of the half-space. 

In the context of the applications already discussed, 
one focuses on a specific resonance and it can be argued 
that a sound figure of merit is afforded by the quantity 
UbT /hiOfiQ. Therefore we compare the different types 
of resonances (for different lengths L) for the same res- 
onant frequency cj^; (fourth and sixth columns). It is 
clear that for specified materials these results for the Q- 
values only depend on the ratios between the beam's di- 
mensions. Thus, for semiconductor heterostructure real- 
izations the corresponding formulas (second and third 
columns) hold all the way from the nanoscale to the 
macroscopic regime and are also applicable to microme- 
chanical resonators. However, while for fabricated flex- 
ural resonators with sub-micron transverse dimensions 
aspect ratios L/d ~ 10 are not uncommorP, for typical 
micromechanical resonators used in MEMS L/d > 100 
rendering the above results an upper bound that would 
be hard to reach — note for comparison that the last 



two entries of the fourth column correspond, respectively, 
to L/X = 49 and L/\n = 24. Finally, we note that for 
given 3D supports the above results are "universal" in 
the specific sense that their only dependence on proper- 
ties of the beam is through /ifc, {r'^)s, and the prefac- 
tor of the TRE dispersion relations quantities that 
specify its low-frequency effective theory. In particular, 
the scaling laws with the beam's length are completely 
general. Furthermore, for a monolithic structure the de- 
pendence on the material's Poisson ratio in the relevant 
range 1/5 < cr < 1/3 is so weak (cf. Sec. IV I that the 
Q-value of a given resonance can be effectively regarded 
as a geometric property. 

The other instance of interest to which Eqs. ([s]), ^ 
apply and that also yields simple scalings, is the case of 
slab 2D supports of thickness dg = X equal to the thick- 



ness of the beam. The results we will obtain [cf. Sec. IV 



Eq. (64l] for the frequency dependence of the Pfj{io) to- 
gether with Eqs. (|2| imply the following 



QnUl ' 



(9) 



which correspond to 1// nois^^^^. We note that the 
finite size of the slab provides a natural infrared cutoff. 
The corresponding Q-vahies (which with the exception 
of horizontal bending scale as the aspect ratio L/w) have 
already been derived in Ref. [33' up to the prefactor C„ 
discussed above. 

Standard fabrication procedures^ normally result in an 
undercut of the support of size djj at least comparable to 
the width of the beam. Our analysis of the abrupt junc- 
tion(s) with the support(s) given in Sec.|lV]and App.[C] 
implies that it is quantitatively correct to use the 3D 
model for the support(s) when djj < d, and the 2D model 
discussed above when djj ^ L. In addition, based on 
heuristic considerations, the 3D results for the Q- values 
given in Table [T] are expected to be qualitatively correct 
for 64 du < L. 

Finally, it should be noted that a substantial part 
of our derivation of the spectral density for an isolated 
resonanc^^ is independent of the geometry. More pre- 
cisely Eqs. ([5]) and ^ can be viewed as a specific instance 
of the following more general relation 



dr^ {u'^ 
s 

X 5[uj 



(10) 



where Uq(f) and (Tq{r) are the displacement and stress 
fields associated with scattering eigenmodes for the whole 
structure, labeled by q [eigenfrequencies i^(<z)], and u^(f), 
cr'j^{f) are the analogous fields for the resonator mode 
(cf. Section|ll]). Here S is the contact surface between the 
resonator and its support (s) and ps is the mass density of 
the latter. As will become clear in Subsec. IVBI the small 
parameter associated with the above approximation is 
\Aj(LUfi)\/2LLij^, where Aj{ujfi)/2 is the support-induced 
shift of the resonant frequency (cf. App. [f|. 
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FIG. 2: Schematic diagram of a microtoroid structure with 
largest diameter D and membrane thickness h (the origin is 
set at the center of the contact area S with the pedestal). 



In all natural scenarios \Aj{ujfj)\/2ujj^, <C 1 arises due 
to the abrupt nature of the junction(s) with the sup- 
port (s) and the condition \/Skpi <C 1 — where kji is the 
typical wavevector associated with the resonator mode. 
The behavior as S* ^ leads to two possibilities: (i) 
the limit is singular, or (ii) it defines a well-behaved res- 
onator geometry that can be described by 3D elasticity. 
The beam geometry falls into case (i) for which is 
specified by clamped boundary conditions at S (cf. Sub- 
sec. II B). In turn in case (ii) u'j^if) should satisfy free 

An exper- 



boundary conditions at the contact surface, 
imentally re levan t example of the latter is afforded by 
microtoroid^^^ (cf. Figure |2| . For such structures an 
heuristic treatment of the pedestal as a beam with adi- 
abatically varying cross section allows us to obtain from 
Eq. ( 10 ) the following 



\/PsEs 

rriRUjR 



R.z 



(11) 



for any axially symmetric isolated resonance. Here the 
resonator mode is normalized so that the normal coordi- 
nate corresponds to the elongation of the toroid's exter- 
nal radius, and mn and Ui?(r) are, respectively, the cor- 
responding effective mass and mode profile (cf. App. [D|. 
Equation ^ then yields for the corresponding Q-value 



(12) 



Typical values for the radial breathing mode of state 
of the art structures ar^^ ojr = 2tt x 50MHz, mu — 
10-"Kg, S = 7r/4x (0.5^)2, and ur^^{0) = 1/2, which 
(for an Si substrate) result in Q w 4 x 10"^. In addition 
to ^/Skfi <C 1, the above approximation assumes that 
there is perfect impedance-match between the pedestal 
and the substrate, and that V~S < h — here h is the 
smallest characteristic dimension of the resonator. It can 
be argued that for pedestals with lengths at least com- 
parable to 277 /kn deviations from this adiabatic scenario 
will only increase the Q- value provided that does not 
coincide with a resonance of the pedestal. The derivation 



II. BASIC MODEL AND OUTLINE 



Normal modes of the supports, the beam, and 
the whole structure 



We will make two simplifying assumptions: (i) the 
supports have typical dimensions that are much larger 
than the beam's length L and for a bridge both of them 
are identical, and {ii) both the supports and the beam 
present refiection symmetries with respect to the xz and 
yz planes (cf. Fig. [ij. The first assumption is adequate 
to describe a wide range of realistic structure^^. The 
second assumption is a theoretical simplification and has 
no impact on our main results [Eqs. ([5]), ([6| and the first 
two columns of Table [l] as will be borne out in App. [C] 
These assumptions allow us to model the supports as 
elastic half-spaces. 

Thus, the whole structure presents reflection symme- 
tries with respect to the x and y axes that are associated, 
respectively, with the operators and Ry acting on the 
space of solutions of the elastic wave equations. This 
allows us to find normal modes of the whole structure 
that are eigenvectors of R^ and Ry. For a given normal 

mode \u) one can generate modes ), \u ^), 

and \u ) with the desired refiection properties by the 

following symmetrization procedure 



\u^i,)^S^i,\u) = -[\u) + p.Rx\u) + vRy\u) + fiiyRxRy\u] 



with fi,i> ~ ±, which guarantees: 



We note that 



E I 

— ± 



(13) 



(14) 



(15) 



Furthermore as in any standard scattering problerrP^, 
we can choose for these modes {\u)} incoming scatter- 
ing states that present simple asymptotic behavior for 
t — > — oo corresponding to the different modes that can 
propagate in the supp ort fs) — i.e. the free modes of an 
elastic half-spac d^^ l '^^ i Thus for each type of mode with 
well defined refiection properties, there will be four types 
of asymptotic behavior, namely: longitudinal bulk waves 
(I), transverse bulk waves polarized along the plane of 
incidence (SV waves), transverse bulk waves polarized 
perpendicular to the plane of incidence (SH waves) , and 
Rayleigh surface waves (s) — also known as SAW. In all 
four cases the corresponding unsymmetrized solutions for 
a free elastic half space can be written in the form: 



(27r) 



[eof 



iq-r 



Aieie 



iqi-r 



AtEte 



^qt■r^ 



of Eqs. (10 1 and (11) is given in Appendix D 



(16) 
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where Qz > [cf. Fig. [T] and 7 = Z, SV, SH, s labels the 
type of mode. The first term corresponds to the incident 
wave and the last two to the reflected longitudinal (l) 
and transverse (t) waves. The dimensionality dj, polar- 
izations e, wavevectors qi/t and amplitudes Ai^t (which 
depend on q and 7) corresponding to each of these types 
of modes are given in Appendix [X] In addition for the 
bridge geometry there will be right- and left-movers and 
therefore eight types of modes. Thus, each mode Wg'^^(f) 
of the support on the left with Qx, Qy, Qz > will gen- 
erate four incoming "right-moving" symmetrized normal 
modes of the whole structure ■Sg^-y_^_iy,R,(f) — with anal- 
ogous relations and definitions for the left-movers and 
in the case of a cantilever. In the following, to simplify 
the notation, we will absorb the discrete indices in the q 
[i.e. Mg,7,;i,i/,R/L('^) Uq{f)] unless otherwise stated. Of 
course, for the bridge geometry we can exploit the reflec- 
tion symmetry with respect to the x-y plane through the 
midpoint of the beam to obtain the left-movers from the 
right-movers so that only the latter need to be calculated. 

On the other hand the beam can be viewed essen- 
tially as a phonon waveguide. Thus it presents a series 
of branches associated with the size quantization of the 
transverse wavevector. An ubiquitous feature of these 
systems is that there are four branches that lack an in- 
frared cut-off: two bending branches (vertical and hor- 
izontal) with quadratic dispersion relations, and a tor- 
sional a nd a compression branch with linear dispersion 
relationd^^EH Each of these branches corresponds to 
a given type of symmetrized modes Iw^jy) — for 

compression (c), ) for vertical bending (u), \u |_) 

for horizontal bending (h) and |u ) for torsion (t). As 

we are interested only in low frequency phonons these 
are the only modes that can propagate inside the beam. 
Henceforth we will identify the discrete indices /i, v of 
the normal modes Uq{r) with the corresponding branch 
inside the beam by introducing a suitable branch index 
P = c,t,v,h (i.e. {HjI'} — > /3). As already discussed, to 
describe these branches at low frequencies one can resort 
to TREP, that consists of an approximation in which the 
small parameter is kd. More precisely it corresponds to 
taking the lowest order terms in an expansion in kd of 
the transverse profile of the mode. The associated dis- 
placements are given in Appendix [B] For a finite length 
beam, in addition to these propagating solutions, it is 
necessary to also consider the role of exponentially de- 
caying evanescent solutions. These are of two types: end 
corrections that decay over a length scale of order d and 
solutions associated with the bending branches with a 
decay length 1/k — the latter can also be treated within 
TRE. 

As we are interested in properties of the beam's motion 
our main task will be to determine Uq{f) inside the beam 
(and a distance d away from the junctions). This prob- 
lem is analogous to the one of finding the electromagnetic 
field inside a Fabry-Perot interferometer. Thus it can 
be approached by considering how the elastic waves that 
can propagate inside the beam are reflected at the abrupt 



junctions (and at the free end for the cantilever) and then 
adding the infinite refiections. In the cases of the bending 
branches for fc < 1/L it is clear that the real exponen- 
tials are not simply end corrections and have to be taken 
into account in this resummation process. Hence there 
are two ways of propagating for the bending modes. This 
reduction of the problem to the refiection at a single junc- 
tion is presented in Section |III[ In addition one needs to 
consider how a wave incident from the support is trans- 
mitted into the beam. These problems of reflection and 
transmission at a single junction are tackled in Sec. |IV[ 
Heuristically, for very long wavelengths the abrupt 
junction is seen from the inside of the beam as a clamped 
boundary condition, and from the outside (i.e. from the 
support) as a free boundary condition. Let us consider 
for example the problem of reflection from the inside. If 
one adopts as an approximation for the displacements in- 
side the beam the solution that corresponds to a clamped 
boundary condition u^{r) it can be proved that the as- 
sociated relative error is at most of order kd (cf. Sub- 
An analogous treatment can be done for the 



IVAl 



problem of transmission into the beam from the out- 
side. It will be proven in Subsec. |IVB| that the associ- 
ated transmission amplitudes can be exactly related by a 
reciprocity argument to the total transmission coefficient 
into the support t^uj) associated with reflection from the 
inside of the beam. This and the aforementioned approx- 
imation by free boundary conditions (cf. Subsec. IVAl 



will allow us to establish simple expressions for the func- 
tions p{lu) and t(lu) valid to lowest order in kd. 

Finally, in Sec. [V] we will relate the environmental 
spectral density /(w) for each resonance with the nor- 
mal modes Uq(f) and establish the fundamental rela- 
tions ([5|, ^ between the functions I{uj) and p{lj), t{lij). 
As will be discussed in the following Subsection, their 
derivation involves the use of a canonical transformation 
that relates a closed quantum system "scattering repre- 
sentation" for the "mechanical Fabry-Perot" (furnished 
by the supported beam) with an open quantum system 
"resonator-bath representation" that singles out one of 
its resonances. 



B. Closed system vs. open system 

To each normal mode Uq{r) with frequency i-u(q) we 
can associate a bosonic annihilation operator b{q). Thus, 
the normal ordered Hamiltonian operator for the whole 
structure considered as a closed quantum mechanical sys- 
tem adopts the simple form 



/ Mg)6t(g)6(g). 

Jq 



(17) 



Here denotes dqxdqydqz and summation over the 
appropriate discrete indices detailed in Subsec. |II A| In 
this diagonal representation the resonances will emerge 
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as poles of the propagator for the displacement field 



u{r) 



2psUj{q) 



Uq{f)b\q) +11.C. 



(18) 



linearly coupled to a fielcP^. The dynamics of the res- 
onator mode is completely determined by the following 
environmental spectral density: 



and in this sense may be regarded as "derived objects" . 
It is important to note that when the beam is character- 
ized by a different density pt ^ ps the above eigenmodes 
of the whole structure are orthonormal in a non-trivial 
metric which defines the scalar product that will be used 
in general (except in Sec. |IV|and Apps. [B]and[C|: 



/(c.)^2^ / \aq)?5[u-u{q)] . 



v*{r) ■ u{r) . 



(19) 



(21) 



This problem is integrable given that H can be diagonal- 
ized by a linear canonical transformation^* that would 
allow to determine the Uq(f) given u^(f), u'^^^f), lob, and 
C(g). We are simply faced with the inverse problem: the 
derivation of Hamiltonian (20 1, with 'S'^(f) and {C,{q)} 



On the other hand in the context of the applications 
already discussed (cf. Sec.|l|, where the system is driven 
out of equilibrium, it is normally convenient to isolate 
explicitly the relevant degrees of freedom of the nanores- 
onator and treat them as an open quantum system. In 
the simplest scenario there is one resonance of interest so 
that the goal is to find a "resonator-bath" representation 
in which the Hamiltonian adopts the form 



satisfying the aforementioned properties, from Hamilto- 
nian (17 1 given the normal modes {uq{r)}. The solu- 
tion is not uniqu^^ but quantities of physical interest 
— e.g. the resonator's frequency dressed by the environ- 
ment, its Q- value and its relative oscillator strength in 
physical spectra — will remain invariant under different 
choices. 

Our analysis of the long wavelength behavior of the 
propagation of elastic waves inside the resonator (cf. Sub- 
and 



sec 



IVAl 



H 



h b'+ b'l 



[C(g)6'((7) + H.c.] 



hL0{q)b'\q)h'{q) 



(20) 



the fact that the resonances of in- 
terest have low frequencies that correspond to wave- 
lengths much larger than d, prompts us to choose for 
M'^(r) an eigenmode of the elastic wave equation satis- 
fying clamped-clamped (clamped-free) boundary condi- 
tions for the bridge (cantilever) geometries. It is straight- 
forward to prove, using the expression of the elastic en- 
ergy density in terms of the displacement field-^-™, that 
for this choice ujn will be given by the corresponding 
eigenvalud^. On the other hand the diagonal represen- 
tation afforded by the normal modes Uq{f) directly yields: 



Lol = / u:\q)\{uq\u'j,)\' 



with off-diagonal couplings C,{q) that are sufficiently weak 
so that the renormalization of the bare frequency lor is 
smaller than the characteristic spacing between the res- 
onances — i.e. the free spectral range of the mechanical 
Fabry-Perot. In such a representation the relevant reso- 
nance will correspond to the degree of freedom described 
by the operators 6'^, b'^ — which annihilate or create a 
quanta in the resonator mode u^(r) — and in this sense 
can be viewed as a "fundamental object" . The discrete 
mode u'j^(r) should be real, localized in the beam, have 
finite averaged elastic energy density per unit amplitude 

given by u;l/2, and satisfy the elastic wave equation ;;nTatao"n oTch'^nce'is'i^'^ 
mside the beam. In turn the b (q) will be annihilation 
operators for a continuum of modes Uq{f) [the environ- 
ment] that have support in the whole structure and are 
labeled after the Ug{r) with which they are in a one to 
one correspondence^. More precisely, the classical state 
Uq(f), evolved with the Hamiltonian density that corre- 
sponds to setting ({q) = in Eq. (20 1, and the classical 



(22) 



The required canonical transformation from the diag- 
onal representation (17 1 to the "resonator-bath" repre- 



biq) 



^(g) _^ I 



state Uq{f), evolved with the elastic wave equation for 
the whole structure, will share the same "free" asymp- 
totic behavior for t — > — oo. The form of the off-diagonal 
term in which the environment only couples to the canon- 
ical coordinate of the resonator Xr reflects the fact that 
the underlying phonon Hamiltonian is diagonal in the 
momenta. This, after integrating out the resulting equa- 
tions of motion for the b'{q) allows to obtain the well 
known Equation^* (| l|). 

The Hamiltonian (|20[) is just an instance of the thor- 
oughly studied problem of a single harmonic oscillator 




{Uq\uR) b'ji 



^R 

Lj{q) 



{Uq\uR)bR 







u{q') ^ 


V ^(9) 


uj{q) 




u{q') 





{Uq\Uq,)b'{q') 



{Uq\u':)b'\q') 



(23) 



and the analogous relation implied for the fo^(q), with the 
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spatial overlap {u^\u'^,) given by: 



Here we have defined, for any smooth function f{uj) 



(24) 



1 9i / 

Afic)^-^ du;'- ^/(c.'), (25) 

and introduced the function 

S{u;)^^lj{u'j,]u,)\'S[u;-u;{q)]. (26) 

We note that this definition remains invariant if we re- 
place the Uq{r) by any other basis of normal modes of the 
whole structure, so that the function S'(w) is a property 
of the resonator mode. In fact, it corresponds for w > 
to the spectral density of the resonator mode's canoni- 
cal coordinate at zero temperature. The choice of ~ie 



in Eq. (24) yields the desired asymptotic condition for 
t —oo, already discussed. The normalization of the 



resonator mode and Eq. ( 22 1 imply two sum rules 



UjS{uj)duj — — and / ui^ S{uj)duj = ■^z^^r , 

2 Jg 2 

(27) 

respectively, which will prove useful below. Here we have 
used that for any smooth function /(w) 



\{u'j,\u,)\'fHq)]^ I dJ-^S{u;)f{^). (28) 



{Uq\u'*,) 



u^,\Uq) can 



and the overlap (24 1 



On the other hand the overlap 
be expressed in terms of 

inserting 1 = ^ij\uq){uq\. In turn, as complex conju- 
gation corresponds to time inversion, is just the 
5-matrix for the corresponding elastic wave classical scat- 
tering problem^^ . The facts that the ^-matrix only mixes 
states with the same frequency and that the resonator 
mode u'j^{r) is real can then be used to establish: 



uj{q') V uj{q) 



{u \u'*,) ^ 



l^{q) / ^(gQ 
uj{q') V Lo{q) 



{u \u'ji){u'j^\u*) 



As[uj{q')] + iS[u:{q')] uj^q) - LJ^q') + 



To proceed, one can first prove that the inverse transformation has analogous form given by 

{u'^\uq,)h{q')+ j ± J^-^/^ {u'^\ul,)b'^{q') 



b'{q) 



, 2 
1 

, 2 



L0{q') 



uj{q) 



{u'q\Uq,)biq') 



1 

, 2 



iu{q') V u;{q) 



{u'q\ul,)b\q'). 



(29) 



(30) 



We have done this by showing that substitution of 
Eq. (|23| and its counterpart for 6^(g) into the RHS of 
Eq. ( |30[ ) yields the identity, which reduces to a straight- 
forward calculation with the help of Eqs. (24l-(29 1. Sub- 



sequently, in an analogous fashion, one can use Eq. (30) 



to prove that the b'{q), b'\q) also satisfy canonical com 
mutation relations. Then, substitution of Eq. (23 1 an d its 
counterpart for b^{q) into the diagonal Hamiltonian ( |l7[ ) 
leads to the desired "resonator-bath" form embodied in 



Eq. (20) with off-diagonal couplings given by: 



aq) 



{UR\Uq) 



2^/ujRUj{q) As[uj{q)] - iS[Lu{q)] 



(31) 



Finally, substitution of this result into the definition (21 ) 



for I{uj) allows us to obtain 



ujR A|(tj) 



(32) 



for Lu > 0. This together with Eq. ( 26 1 reduces the deriva- 



tion of the environmental spectral density to the calcula- 
tion of the overlaps {u'j^\uq) and implies that it is entirely 
determined by the choice for the resonator mode. We 
note that both the Uq{f) that diagonalize the Hamilto- 
nian (20 1 and the Uq{f) that solve the inverse problem for 
a given u'fj{f) are fixed once their asymptotic behavior is 
specified. 

It is worth noting that within the above exact treat- 
ment the {Cq} describe the coupling to an environment 
that — barring selection rules arising from the symme- 
tries — includes the other high-Q resonances of the 
beam coupled to the supports. Thus the naive expecta- 
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tion that /(w) be smooth on the frequency scale lob. will 
fail in a neighborhood oi uj = ujn ^ 'jJr where the envi- 
ronment is structured and /(w) may exhibit interference 
effects. This issue and the natural ultraviolet cutoff for 
/(w), which is set by the transverse dimension d, will be 
discussed further in Sec. |V] and App. |F] 

III. REDUCTION TO A SINGLE JUNCTION 

It is clear from the analysis in the previous section that 
we will eventually need an expression for the overlaps 
{u'jj\uq). As u^(r) only has support inside the beam, 
knowledge of the normal mode Uq{r) is only needed in 
that region. Furthermore, it is straightforward to realize 
that the end corrections (cf. App. |b]) will only yield con- 
tributions to {u'f^\uq) that are higher order in d/L so that 
we may focus on the four branches that can propagate 
at low frequencies. As already discussed in Subsec. II A 



given our 



our system can then be regarded as a Fabry-Perot in- 
terferometer for elastic waves and for the bridge we may 
concentrate on right-movers. It proves useful to explic- 
itly separate the transmission amplitude tq for propa- 
gation inside the beam that would ensue for the corre- 
sponding problem of a single junction (cf. Sec. IV I; since, 
given our model of the supports (cf. Subsec. II A I, in the 
limit d/L — > the combination Uq(f)/tq for f inside the 
beam only depends on q via ujq and the branch index f3. 
For each branch the latter frequency will correspond to a 
wavevector k(q) for propagation inside the bearrP^I and, 
as described in App.|B] the mode profile Uq{f)/tq can be 
expressed in terms of an effective one-dimensional field 
4>k{q),i3{z) [cf. Eq. B8 — henceforth in this Section we 
will omit the g-dependence and the index (3. 

In the cases of compression and torsion only traveling- 
wave solutions intervene (cf. App. [B|. Therefore 



0fe(z) = A+(fc)e^'=- + A_(fc)e 



-ikz 



(33) 



where we take < z < L, and A±{k) are the amplitudes 
of the right and left moving components of (j)k (cf . Fig. [ij . 
In order to determine the above amplitudes we use the 
usual procedure for Fabry-Perot interference, i.e. we de- 
rive the amplitudes from the infinite sum of contribu- 
tions generated by consecutive reflections at the ends. 
We find that they obey the following linear system of 
self-consistent equations: 



(34) 



where we have defined A^^^^ as the amplitudes for reflec- 
tion at the right /left ends of the beam, and the flrst entry 
on the RHS corresponds to the displacements generated 
by the incoming wave from the support. This system of 
equations is valid for both the bridge and the cantilever. 
For the bridge geometry we can use the reflection sym- 
metry with respect to the x-y plane through the mid- 
point of the beam to prove that the reflection amplitudes 




are related by A^'^ — 7.(^)5*2'=^', Furthermore, 
choice of origin, dimensional analysis implies that for a 
speciflc support material and type of beam the reflection 
amplitude is just a function of kd (analogous consid- 
erations will apply below to , with S,r] = A,B, and 
to the ratio b). 

On the other hand, as discussed in Sec. Ill A| and 



App. |BJ in the cases of the two bending branches one 
needs to also consider evanescent solutions that decay 
over a length scale of order 1/k. In order to include 
these in our analysis, we generalize the above procedure. 
The functions (fik corresponding to the bending modes 
will have four contributions with respective amplitudes: 
A±{k) for propagation to the right/left, and B±{k) for 
decay to the left/right. Hence we write 

c^k{z) = A+(fc)e*'^"H-A_(fc)e-'^'"-fB+(fc)e'^'"+B_(fc)e^'=" 

(35) 

In formal terms the resummation of the successive "re- 
flections" can be viewed as an iterative procedure in 
which the nth contribution when added to the n — 1th 
matches the boundary conditions at the right/left ends 
for odd/even n (here n — 1,2... with n — Q correspond- 
ing to the solution for a single junction). Thus, "reflec- 
tion" at one end of the beam will generate both types 
of contributions, i.e. there are now eight reflection am- 
plitudes. We can define, for example, as the am- 
plitude for a propagating mode to be reflected at the 
left junction into a decaying mode and similarly for the 
other amplitudes. Hence the system of self-consistent 
equations obeyed by the amplitudes which specify (l)k{z) 
generalizes to 



/-I 





' AA 


^BA ' 









-1 


JL) 
' AB 


JL) 
' BB 


B_ 


-b 


' AA 


' BA 


-1 










^^AB 


' BB 





-I J 


\bJ 


Vo/ 



(36) 



where b corresponds to the decaying contribution gen- 
erated by the incoming wave from the supports. More 
precisely it is defined as the ratio of the decaying con- 
tribution's amplitude to tq for the corresponding prob- 
lem of transmission from the support into the beam at 
a single junction (see Sec. |IV[ ). For the bridge geom- 
etry we can reduce the eight reflection amplitudes to 
four using again the corresponding reflection symme- 
try. 



This yields 
"^''^ - r(i']e-(i-^)'=^, and r 



JL) ,2kL JR) 
' AA^ ' ' AB 

(R) - SL) 2kL 
BB ^ ' BB^ 



' AB^ ' 



' BA~ ' BA"" 

Finally, we note that if one allows the ratio b to depend 
on g, the results of this Section for the modes Uq{f) inside 
the beam (and a distance ^ d away from the junctions) 
are applicable to the extent that contributions of order 
exp[— L/(i] arising from the end corrections are negligible 
so that Eqs. (33l-(36l are completely general for k{q)d ^ 
1. 
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IV. TRANSMISSION THROUGH A 3D-1D 
JUNCTION 



A. Approximation by clamped and free boundary 
conditions 



We now turn to the analysis of the small kd behav- 
ior of the transmission amplitudes tq, btq and reflection 

amplitudes r^^j, i"lf/hSri (^ith 6,1] — A,B) that charac- 
terize a single 3D-1D abrupt elastic junction (cf. Sub- 
sec. II A I. To this effect we consider incoming eigenstates 
that present simple asymptotic behavior for t —oo. In 
the case when the incoming wave is incident from the sup- 
port these eigenstates \uq) will correspond to the different 

free modes of an elastic half-space u''q^\f) already dis- 
cussed in Section |II A[ Alternatively, for incidence from 
the beam, they will correspond to the left-movers of an 
infinite beam. Once again we consider modes with well 
defined reflection symmetries. This implies for the eigen- 
modes incident from the support the same labels g, 7, /? 
as for the modes of the whole structure while the modes 
incident from the beam are specified by k, [3. 

We focus on a specific eigenmode u{f) and for simplic- 
ity omit its labels. Inside the beam [i.e. for z > and 
(x, y) in a beam's cross section] we decompose u{f) into 



u{r) = Au^{r) + u*(r) 



(37) 



where (f) is an approximation to the displacement field 
specified by taking at z = 0: (i) clamped boundary con- 
ditions when the incoming wave incides from the beam 
or (a) displacements specified by the corresponding so- 
lution for the free elastic half-space [i.e. u*{x,y,0) = 
S u^^\x,y,0)] for incidence from the support [cf. Eq. (13) 



and App. |A]. The problem of finding u{r) inside the 
beam can then be formulated as an integral equation for 
Au+(r): 

Au+{f)^ [ dr'^G{r~f',uj)-F-[Au+{f') + u^{r')], 
Js 

(38) 

where S is the beam's cross section at the origin, r — 
{x, y, 0) e S, and z' > 0. Here F and G are second rank 
tensors. The former is given by the linear differential 
operator that maps u{f') onto z ■ cr{f'), where (T{f') is 
the induced stress tensor, and the latter by the retarded 
Green's function of the free elastic half-space harmoni- 
cally forced at its boundarj^H. More precisely, Gij{f,uj) 
is defined as the i-th component of the outgoing dis- 
placement field generated at point f by the harmonic 
stress source [frequency w(fc)] with amplitude specified 
at z' = by: (T,,(a;', y', 0) = 6{x')S{y'), ai,(x\y',0) = 
for I ^ j. In order to establish that Eq. (38 1 specifies 



the solution for z > we just need to show that the ex- 
pression ( 37 ) satisfies continuity of the displacement and 



the stress with a solution in the half-space {z < 0) hav- 
ing the appropriate asymptotic boundary conditions. We 
define Au_(f) as the extension of the RHS of Eq. (38) 



for z < and arbitrary x,y. In case (i) this function di- 
rectly gives the required solution in the support, while in 
case (m) the latter is afforded by Su^°\f) + A-u_(f). In 
both cases the continuity of the stress follows trivially by 
construction whereas Eq. ( 38 1 enforces the continuity of 
the displacement. It is understood that both Au^{f) and 
(f) are linear superpositions of the low frequency har- 
monic solutions |t;i™''(fc)) of a semi-infinite elastic beam 



(cf. App. |B| with frequency Lo{k) (specified by k(q) for 
incidence from the support). Thus u{f) satisfies the elas- 
tic wave equation for z > and the only traveling wave 
contributing to A?2+(f) is the right-mover corresponding 
to k, (3. This yields 

Aw+,/3(f, k,d) = Y^ c^m(fc, d) 

rn 

X Afi^{x, y, d, K/3™)e-'"''"('=''')^ , 
u*,f3if, k, d) = Uin,i3{f, fc, d) -f ^ ^) 

m 

X A^,^(x,y,d,K0„)e-^'""('=-'*)^ (39) 

Here we have reintroduced the "symmetry index" /3 and 
eliminated the frequency in favor of k. The ampli- 
tudes in the different harmonic solutions are given by 
cp,m{k,d) and c^*^(/e,d); while Af}^rn[x,y,d, K,(}^rn{k,d)] 
and Kf3^m{k, d) are, respectively, the corresponding trans- 
verse profiles and complex wavevectors of the lat- 
ter. The displacement field incident from the beam 
Uin,0{^, k, d) vanishes for case (ii) whereas for case (i) it is 
given by Ai^fi[x,y,d,ik]e~'^''^ for propagating modes and 



A. 



v/h,l 



[x,y,dr 



j^]gK„/fi,iz £qj. large decay length 



exponentials associated with the bending branches {de- 
caying modes). This "incident displacement" and the 
terms with m = and with (3 = v/h,m = 1 yield the 
TRE part of the solution, while j3 = c,t with to > and 
f3 = v,h with TO > 1 correspond to the end corrections. 
These are characterized by 5ft[K^^m(A;, d)] ^ 1/rf whereas 

Kfifiik, d) = -ik Hv/h,i(k, d)^k{\ + 0[kd]) . (40) 

To extract the small kd behavior of the transmission 
and refiection amplitudes we first prove that 



|A^+(r)| ^ 
\u{f)\ - 



0[kd] for z>d. 



(41) 



For this analysis we eliminate k and z in favor of kd 
and kz, respectively, which are then treated as indepen- 
dent variables (henceforth we omit the resulting d depen- 
dence of non-dimcnsionless amplitudes). If one substi- 
tutes Eq. (pTl) into |Au4.(r)|/|u(f)| one can deduce that 



|Ai2+(r)| 
|^.(r)| 



< 



0[kd] 



|Aa+(f)| 
\uir)\ 



< 



0[kd], (42) 



so that it suffices to prove the LHS of Eq. (42 1. Further 



more, for z ^ d the contributions of the end corrections 
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are exponentially suppressed so that it suffices to ana- 
lyze the TRE amplitudes [cf. Eq. (39l]. In fact it can be 
proved that the latter satisfy 



,(*) 



(fed) 



\Cf3,m{kd)\ 



< 



< 



{0[kd] for case (ii) and P = t 
0[1] otherwise 

io[{kd)^] for case (m) and /3 = t 
[©[/cd] otherwise 



(43) 



(44) 



which then directly imply the LHS of Eq. ( 42 1 



First we establish the behavior (43 1 of the starred am- 



plitudes which amounts to a rigorous derivation of the 
recipes used in TRE to specify the boundary conditions 
for the effective one dimensional field (p{z). To this effect 
we substitute Eq. ( 39 1 into 



u^^f3{x,y,0,k,d) 



for case (z) 

S i3ulj'\x , y , 0) for case (ii) 



(45) 



and take on both sides the spatial averages {■ ■ ■)s ^^^d 
(• ■ ■)ang,s that correspond, respectively, to the displace- 
ment of the center of mass and the spatially averaged 
angle 9 [as defined in App. |b] Eqs. ( |b11| ), ( |B13| )] for 
the cross section S (henceforth we omit the latter label). 
If we consider each component separately we have six 
equations that are linear in the amplitudes c^*^j(fcc?) with 
inhomogeneous terms arising from -~Uin,p{x,y,0,k^d) in 

case (i) and S i3u'^\x , y , 0) in case (ii). The reflection 
symmetries imply that in the cases of compression and 
torsion, respectively, only the equation corresponding to 
(uz) and to 9z does not vanish trivially, while for vertical 
(horizontal) bending the same applies to the two equa- 
tions provided by (ux) {{uy)) and 9y {9x)- In each of these 
equations we solve for the TRE amplitudes in terms of 
the end corrections and the inhomogeneous terms. The 
small kd behavior of the resulting expressions can be 
extracted using the following properties of the modes' 
transverse profiles: 



for m — Q 
otherwise 

for m = 0, 1 
otherwise 

for m = 
otherwise 

for TO = 0, 1 
otherwise 



( / d/'G (f - f', u[k]) ■ F- [x', y', d, Kf,,,^{k, d)] 

Js 




-1^13, mik 



0[{kd)Pi'] for /?, TOG TRE 
O [{kd)^Pf] otherwise 



( / d/'G (r - f', iu[k]) ■ F {A^,™ [x', y', d, Kp,^{k, d)] 



xe 



0[{kd)P''] for /?, TOG TRE 
O [{kdfPf] otherwise 

(46) 



which imply that the contributions of the end correc- 
tions scale at most as the inertia (i.e. as cj^). The above 
Eqs. (46 1 follow directly from the universal properties 



of the end corrections discussed in App. |Bj (Subsection 
B 2 1 and the small kd behavior of the TRE solutiond^ 
(cf. Subsection |B 1[ ). Thus, with the help of the Taylor 
expansion of Sf}^\x,y,0) at the origin, Eq. (13 1, and 



Eq. (B12l we obtain for compression and torsion: 



.(*) 



— 1 for case (i) 

2\/2^ui°^(0) for case (a) 



J*) 



^{kd) 




z • V X 4°^(0) + 27,ui°J,,(0)l {ii) 

(47) 



(0) 



where 7z = {ly — Ix) I Iz and Uij is the strain tensor. In 
the case of the bending branches for case (z) we get: 



4}ko(0) 



i for propagating u 

— (1 — i) for decaying u 



Uin,/3 



4}li(0) = 



-( 

— z 



1 + z) for propagating ui^^p 
for decaying Uin,/3 



(48) 



while for case (zz) we obtain: 



1 1 
±i Tl 



.(*) 







(49) 



where the first equation corresponds to the average dis- 
placement (ux/y) and the second one to the angle 9yfx- 
The fact that the latter does not contribute to lowest 
order in this case is a consequence of the linear versus 
quadratic dispersion relations that characterize the prop- 
agation of the relevant modes (/? = u, h) in the support 
and the beam, respectively. Equation (49 1 yields 



(*) 



o(0) = ^/2^(l-z)4%^^(0) 



v/h 



i(0) 



^(*) 

'^v/h,0 



(0) 



^i. (50) 



An analogous procedure can be followed to derive 
Eq. (44 1. We substitute instead the decompositions 
(39 1 into Eq. (38) and now the averages (•■•)s and 



(• • •)ang,s yield linear equations for the un-starred ampli- 
tudes with inhomogeneous terms arising from the starred 
ones. Then Eq. (44 1 follows from Eqs. (46l-(50l com- 



pleting our derivation of Eq. (41). Furnished with the 
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latter it is clear that the lowest order contributions in 
kd to the transmission amplitudes tq{kd), b{kd)tq{kd) 
and reflection amplitudes r^^l{kd), r^^]^^^(fcd) can just 

be extracted from M*(r). Thus from Eq. (|47[ ) we ob- 
tain for both compression and torsion r^^^(O) = — 1. 
In turn, Eqs. (|48|) and (|50| yield for the two bending 



branches: 5(0) 



(0) 



(o) = -(i-*) 

,(0) 



and r^Bsi^) = ^* 



AS(0) = -(1 + 

The lowest order 



contribution tq to each of the transmission amplitudes 
tq is provided by the corresponding approximation for the 



starred amplitude c^^l^kd) [cf. Eqs. (|47| and (|50|]. The 



needed values at the origin of u'^^ (r) and its derivatives 
are straightforward to obtain from the expressions given 
in App. [A] We note that the symmetries of the half-space 
imply ufl^,^^^^{0) = ufl g^^^^_^ jO). 

Finally, for the cantilever geometry a procedure analo- 
gous to the above derivation of the starred amplitudes for 
case (z) yields the standard TRE recipes for the reflection 
at the free end. The latter given our choice of coordinate 



origin lead to: 
torsion, and 



„(i?) 



AA 



(0) = 



i2kL 



for both compression and 



' BA 

bending branches. 



(0) 



(0) = (1 - ^)e(^-l)^■^ r^^l{Q) = -^e 



pi2kL 



^i1](0) = (l + *)e(-i) 



DfcL 



-2kL 



for the two 



B. Reciprocity relations 

Reciprocity is simply the relation = 
for a standard scattering problem in a time reversal in- 
variant theory^^^^. Here v and u €U label freely 
propagating asymptotic states belonging to the sets of 
interest V and C/, "— " and denote respectively in- 
coming and outgoing scattering states and denotes 
the time reversal operation. The inner product "(|)" is 
assumed to be preserved by the time evolution - i.e. the 
underlying theory is "unitary" . However this preserved 
inner product need not necessarily be the usual overlap 
as in quantum mechanics where the "unitarity" corre- 
sponds to the preservation of probability. In particular 
in our case of elasticity theory it is defined so that {■w\w) 
corresponds to the energy carried by the solution w and 
the unitarity corresponds to energy conservatiorPSl This 
relation is quite general but proves to be especially power- 
ful when F ^ [/, in the sense that aside from possible dis- 
crete indices, the available free eigenstates are essentially 
equivalent. In our specific context an example satisfying 
this last requirement would be two different rods (cor- 
responding to V and U) joined at an abrupt junctioiP^l. 
Reciprocity directly implies that for unit incident power 
in a traveling wave € U (with wave vector kjj) the power 
transmitted into a traveling wave € V (with wave vec- 
tor kv) is equal to the power transmitted into the wave 
corresponding to —ku for unit incident power in the one 
characterized by —ky. 

However our model for one junction does not satisfy 



the above "asymptotic equivalence" since it involves cou- 
pling a ID system (the beam) to 3D and 2D continua cor- 
responding, respectively, to the bulk and surface states 
(SAW) of the support. To overcome this difficulty we first 
consider a model of the junction for which: {i) the sup- 
port is characterized by some finite dimension D, (ii) the 
support states under scrutiny are equivalent to a phonon 
waveguide so that the "asymptotic equivalence" is sat- 
isfied, and {in) for D oo the support tends to a free 
elastic half-space. Then we apply the reciprocity rela- 
tions for finite D and finally take the limit D —f oo. 
For bulk states (i.e. q with 7 £ {Z,SV,SH}) a suitable 
"finite support" is afforded by another beam of square 
cross-section (side D) subject to periodic boundary con- 
ditions on the external faces. On the other hand for 
surface states (i.e. 7 = s) a suitable construction is given 
by a slab of thickness D subject to periodic boundary 
conditions at the external semi-infinite horizontal faces 
and a free boundary condition at the finite vertical face. 
In both cases the decoupled support is exactly solvable 
and the support states of interest are equivalent to a 
phonon waveguide whose branches we index with a sin- 
gle label j. Thus the free support eigenmode^^l read 
\u'j^\k',D)) with dispersion relations u!j{k',D) where k' 
is the wavevector along the waveguide's axis. If we define 
the amplitudes tp_j{k' ,D) such that 



ufl{k\D,f) 



tp.,{k\D)vf[kp,j{k\D),f] (51) 



asymptotically for z —> 00, where \v^^\k)) are the freely 
propagating TRE beam modes (cf. App. |B]) with fc > 
and kpj{k',D) the wavevector in the beam's branch (3 
that corresponds to the frequency Lo.j{k' , D); we have 



lim \ufl{k',D))^\uq) 



lim ti3,j{k', D) = tq . 



(52) 

Taking into account that the modes we consider are nor- 
malized (in the standard Euclidean metric) in all ID, 2D 
and 3D cases, it is then straightforward to realize that 

lim [dk'Y^ |i^,,(fc', D)\^6 [u; - u,{k\ D)] 

= [ \tq\^S{u;-u;q) (53) 



where the {q} run only over modes with the correspond- 
ing P and the integration includes summation over the 
appropriate discrete indices. On the other hand reci- 
procity directly implies 

(54) 

Here we have introduced Tp\k,D) as the transmission 
coefficient into the support branch j for a traveling wave 
of type p incident from the beam with wavevector — fc and 
Pp{iu) [^^^(i^, D)] as the power carried by the normalized 
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free waveguide mode with frequency uj and branch index 
/3 The functions Pf3{uj), Pj{uj,D) have a universal 
expression in terms of the mass density of the respective 
waveguide and the corresponding dispersion relation. On 
the one hand, in complete analogy to the equivalent sce- 
nario for the electromagnetic field, the power carried by 
the waveguide mode is given by the product of the group 
velocity and the corresponding energy per unit length 
(averaged over a period). On the other hand the latter 
can be expressed in terms of the mass density and the 
frequency by first considering the harmonic theory of the 
underlying microscopic discrete lattice for a finite waveg- 
uide of length 27r/fc subject to periodic boundary condi- 
tions at its ends and then taking the continuum limit. 
Thus we obtain 



correct to lowest order in the reciprocal of the aspect ra- 
tio d/L only depend on the amplitudes through these 
quantities pp^uj). On the other hand it is clear that if we 

substitute into Eq. ( 56 1 the approximations t'"q^ given in 
Subsec. IV A the pp{uj) determine the transmission co- 
efficients r^[fc^(a;)] correct to lowest order in kd. It is 
important to highlight that ppiuj) is amenable to reduc- 
tion to a property of the free support at the origin that 
directly relates to its DOS or, in the case of torsion, to 
its vacuum spectrum for the angle [as defined in App. |Bj 
cf. Eqs. (Blll-(B13l] and geometrical properties of the 
free beam (S*, J^, and 7^). In fact we can consider the 
displacement field operator for the free support v!'^\f) 
and decompose it in terms of its normal moded^ Mg"-' (f) . 



4^^"" dfc 



1 



2 



Then, the latter decomposition and Eqs. (47 1, (50 1, and 



4^^^" dfc 



(55) 



(561 yield 



where pb {ps) is the density of the beam (support). We p^^uj) = ^^S* / di e*'^*(4"Ho, t)4°)(0, 0)) 



can substitute Eqs. ( 54 ) and ( 55 1 into the LHS of Eq. ( 53 ) 



and perform the integration, the summation and the limit 
to obtain for each branch (3 the following reciprocity re- 
lation: 



—00 
00 



p,(c.) = / dte^- lim(^(o)(t)^(")(0)) 



p„/,(c.)=^5/ dte-'(^^7^(0,i)^^7JO,0)), (57) 



Pb 

Ps 



(56) 



where Tp(k) is the total transmission coefficient into the 
support for a traveling wave of type /3 incident from the 
beam with wavevector —k, and the {q} run only over 
modes with the corresponding f3. 

Finally we note that the above derivation can be ex- 
tended to the cases considered in App. [C] by suitable 
modifications of requirement (iii) so that in all cases 
the limit D ^ 00 yields the support under consideration 
("3D asymmetric" or slab). Naturally, in the asymmet- 
ric cases the labels q will no longer relate to the beam 



where we use Heisenberg operators and (. . .) denotes 
the vacuum expectation valu^^. This connection abets 
the interpretation of Ppiuj) as an effective environmen- 
tal DOS set forth in the Introduction (cf. Subsec. I A I. 
We note that the latter quantity has dimensions of lin- 
ear density of states and that for compression and bend- 
ing it is proportional to the DOS of the support times 
the area of the beam's cross section. We separate the 
contributions of each type of support modes 7 so that 
Pp{lo) = Pi3,'y{^)- In turn, given the symmetrization 



branch index (3 and in Eq. ( 56 1 we will have t, 



t. 



with the running over all the support modes. In the 
case of a SWNT for which we use the shell "continuum" 
model Pb should be replaced by the surface density of 
graphene ac, while for a nanowire or SWNT for which 



of the modes, the wavevector integration in Eq. (56) only 
involves qx,Qy,Qz > while the symmetries of the half- 
space imply for the two bending branches {t^^l ^ ^ fi\ — 
where we use spherical coordinates for 



\t 



(0) 



,(0) 



g,^,7r/2 — (^,7,^) I 

the wavevector q as in App.[^ One can then perform the 



the underlying model for the 
equate ansatz reads: pb Ub/Nc (cf. App. |C|. In the 
case of a thin plate geometrjISS both ps and pb should be 
replaced by the surface density. In fact it can be argued 
that Eq. ( 56 1 does not depend on any specific proper- 



(fc)) is discrete the ad- substitutions uj = Cjq,v = cos9, integrate over uj and 



ties of the junction or the support and only relies on the 
phonon transport being ballistic. 



Transmission coefficients for each branch 



and eliminate c-v in favor of Ct to obtain: 



PcfiuJj = 3- Mc,7(a)w 

Ps 



Ps 4 



-(A) , 2~{S)f s 



. , Pb'^S _ I \ 2 

Pv/h,-/[t^) = 3- w«//i,7(a)^ • 

Ps Ct 



(58) 



We turn now to the evaluation of the leading contri- 
bution in kd to the LHS of Eq. ( 56 ) for each branch [3 



which we define as pp(uj). On the one hand we will find 
in Subsec. VB that the force spectral densities /„^^(w) 



Here we have introduced the following dimensionless 
constants and functions of the ratio a = (ct/ci)''^ = 
(1 — 2o-s)/2(l — CTs) for the supports' material {as is the 
corresponding Poisson ratio) : 
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Uc,i{a) = -: — / dv 



2n 



(1 -2a + 2af2)2t,2 



[4a3/2\/l ~a + av^{l - v^)v + (1 - 2a + 2av^y' 



Uc,sv(a) 



dv 



(1 - a - w2)(l - w2)w2 
16(1 -a- i;2)(l - u2)„2 (2^2 - 1)4 
(a- l + w2)(l - t;2)v2 



[4Va- l + t;2(l - u2)^; + (2i;2 _ i)5 



'"c,SH=0 Uc,s(a) 



«l3v(") 



"t,SH 

M-«/ft,sv(a) 



ma) 



r- ^ 1 - aa)/2 



v7/2 /-l 



du ■ 



(1 - a + au2)(l - 1,2)2^2 



47r ./o [4a3/2Vl - a + aw2(l - v'^)v + (1 - 2a + 2aw2)2] 

1 f^^, (l-w2)2„2(2^;2_i)2 

2^2 - 1)4 

n 2 



167r Jo 16(1 -a- w2)(l - t;2)2„2 (2^2 - 1)4 



1 

24^ 
1 



1 

16^ 



dv 



(1- z;2)z;(2w2- 1) 



4Va- 1 + 1^2(1 - i'^)?; + (2^2 - 1)2 



dv 



4a5/2(l_Q, + Q,„2)(^_„2)^2 



47r7o [4a3/2Vl - a + aw2(i _ ^;2)„ + (1 - 2a + 2aw2)2 



47r 



/l-Q 



dw 



(2t;2 - l)2i,2 



16(1 - a - W2)(l - w2)2i;2 + (2i,2 _ 1)4 



47r 



dv 



(2z;2 - l)2i,2 



1=^ [4Va - 1 + w2(l - w2)„ _|_ (2^2 - 1)2 



47r 



Uv/h,s{a) 



(59) 
(60) 



(61) 



(62) 



where ^(a) is the ratio of the velocity of propagation 
for surface waves to Ct and is a function of a that is al- 
ways less than unitjEsI "w^e note that for compression and 
bending the tq have a non-vanishing limit for A: — > so 
that the frequency dependence of the quantities p follows 
directly from the density of states of the 3D supportP^. 
Whilst in the case of torsion there is an extra factor of 
ti;2 given that the t^q^ scale as the derivatives of the dis- 
placement field . In all cases the contribution of S V 
support modes has two distinct terms: one corresponding 
to polar angles 6 below the critical angle arccos VI — ct 
for which there is a reflected longitudinal wave, and an- 
other corresponding to angles for which the longitudi- 
nal component is evanescent. The corresponding results 
for the effective environmental DOS for nanotubes and 
nanowires can be obtained from Eq. (|58[ ) via the ansatz 
PbS pb, Pbiz lJ'b{'r'^)s (cf- App. IcF Finally, we can 



use Eqs. (56 1, (58 1- (62) and the TRE dispersion relations 
^ to obtain after summing over the index 7 the follow- 
ing expressions for the transmission coefficients into the 
support 



X 1/2 



/3 _ 

Tv/h{k) = 4u„/,j(a) ( y— 

Tt{k)=A 



1-a E, 
3 — 4a Ebly/x 



Pb 
3/2 



a E. 



1 

12^ 



'y-'^^ #'(a) 



PbS 



3-Aa C 



3/2 



PbIz 



5/2 



(63) 



where we have also substituted the definition of 7^ 
used Ct — \/ Es/2ps{l + ag), and introduced up (a) = 



16 



0.15 - 




FIG. 3: Dimensionless displacements Uc (dashed line), Uy/h 
(solid line), and dimensionless angle u[^^ (dotted line) as a 
function of the ratio a = {ct/ciY . For typical materials of 



interest a a ~ 1/3 which yields Wc(l/4) 
0.019, and u„/h(l/4) = 0.12. 



0.13, u^'''(l/4) 



'J2'y^P,-yi'^)- The latter functions are plotted in Fig. js] 
for all ratios a corresponding to physically allowed posi- 
tive values of the Poisson ratio cTj . 

Finally, we consider the frequency dependencies and 
the scalings with the transverse dimensions for the anal- 
ogous problem of an abrupt junction between a rectan- 
gular beam of width w and a slab (support) of the same 
thickness t ^ w = d [cf. App. |C]. The analogous expres- 
sions for the starred amplitudes yield for the correspond- 
ing effective environmental DOS 



w 

t 



WW 



(64) 



which via Eq. ( 56 ) allow us to recover the following scal- 



ings with kd for the transmission coefficients: 
Ti3^h{k) ^ wk Th{k) ^ (wfc)^ 



(65) 



already derived in Ref. [33] by an alternative method 
(cf. the next Subsection) . Here we have used that {r'^)s ^ 
and C/Iz ^ t^/w^ for t < w. 



D. Multipole expansion method 

The method followed above to derive the transmission 
coefficients relies on using the transmission amplitudes 
from the support into the beam — which given the results 
of Subsec. HV Al can be extracted from the free modes of 
the supports — and then exploiting the reciprocity re- 
lation (56 1. The naive expectation is that an alterna- 



tive method should be afforded by using the far field of 



Au_(r) [cf. Eq. (38l] corresponding to incidence from the 



beam [case (i) in Subsec. IV A to determine the power 
across the surface of a hemisphere with r oo per unit 
incident power. 

The displacement field A?2_(f) can also be viewed as 
the retarded solution for the problem of the free elas- 
tic half-space harmonically forced at its boundary by 



the stress source (Jiz{r) = i - F ■ [Au+{f) -|-u*(f)]. As 
this radiation problem involves a source with a typical 
dimension d much smaller than the wavelength associ- 
ated with its frequency it is natural to use for Au^(f) 
its multipole expansion, whose successive moments will 
yield terms of increasing order in kd. More specifically, if 
we consider the asymptotic behavior of the correspond- 
ing Green's function^^ Gij(f, aj)for r — > oo and the 
nth moment of the stress source^ 



(n) 



(66) 



with ijn — 1,2,3, n ~ 0,1,..., one can show that 
the corresponding contribution to the displacement field 
Au'f^\f) scales as 



-in) 



(f) 



q]^{k)M^"'> 
£Lr 



g*?- > 1 



(67) 



where we have introduced the norm 
{E{^}[<l...i„+jV(« + I)}'/'- Naturally the source 
aiz (r) should also be expanded in powers of kd. However 
our treatment of the solution inside the beam (cf. Sub- 
sec. IV A[ ) only provides us with the contributions arising 
from the TRE part of u*(r) whose order is lower than 
the inertia [{kd)'^^'^]. Thus the multipole expansion will 
only be useful if it can be established that the former re- 
sult in contributions to the transmission coefficient that 
dominate over those corresponding to the end corrections 
and to higher orders of the TRE branches. The relative 
orders follow directly from considering the expansion in 
kd of the RHS of Eq. ( 67 1 . A straightforward analysis 
implies that in general (i.e. without invoking any sym- 
metries) this is non-trivial du^^ to the role played by 
the end corrections' symmetric part of the restriction of 
M^^l^ to ii, Z2 = 1, 2 [M('^)] and by their moments M^^). 
We note that the leading contribution to the far field 
of the compression, torsion, and bending TRE branches 
scales, respectively, as k, k'^, and k^. In the cases of com- 
pression and bending these correspond to point sources 
given, respectively, by a normal force and a tangential 
force applied at the origin. On the other hand for tor- 
sion the leading contribution corresponds to a normal 
torque applied at the origin only when the corresponding 
M^'^) vanishes. For the bending branches the reflection 
symmetries are enough to cancel M^'^^ and <7"(fc) ~ k^ 
implies that higher order moments are irrelevant. In 
the case of the torsional branch it is sufficient to have 
M(^) = M(2) = which can be guaranteed by also re- 
quiring symmetry under the swap x y (i.e. if the cross 
section is symmetric under rotations by tt/2 around z). 
In contrast for compression augmenting the symmetry to 
that of a cylinder does not guarantee the cancellation of 
the diagonal part of M^'^^ which could potentially give 
contributions to the far field of order k. Thus, whilst for 
bending and torsion the multipole expansion approach 
will provide a non-trivial corroboration of the exact form 



of Eqs. (63 1 — and together with Eqs. (58 1 of the red 
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procity relation ( 56 1 — for cross sections with the appro- 
priate symmetries, for compression it in principle only 
allows to check the fc-dependence of these relations. 

Energy conservation implies that the evaluation of the 
contribution to r of the aforementioned leading multi- 
pole of the TRE branch can be done at the origin. In 
general if we consider a harmonic point source given by 
the superposition of a force -Fout and a torque Mout ap- 
plied at the origin we have for the power radiated into 
the half-space averaged over a period 



f:. 



Sout(0)-HM* 



t(0)} (68) 



where 9ont{0,t) = lims^o ("out (j^, i))ang,s-- with S a 
circle on the free face at z = [cf. Eq. ( |B12| ] — and 
Wout(^ji) is the generated displacement field. If we now 
express the lowest order in kd of the incident power as- 
sociated with Uin{r,t) [cf. Subsec. IV A in terms of the 
corresponding total force Fin and total torque Min using 
the TRE transverse profiles given in App.lBlwe get 



[f;,, ■ Sout(O) -t- M*„t • lut(O)} 



F. 



i5R{F;„.^i„(0) + M*„.L(0)} 
Our results for the refiection coefficients I'c^tiO) 



(L) 



(0) imply 



Fout 




Fin 





P = v,h 
P = c 



Mout 









P — v,h 
4 f3^t 

(70) 

the properties of the 
TRE solutions and the definition of the Green's function 
G it is straightforward to establish for the leading order 
of the transmission coefficients 



From the above Eqs. (69 1, (70 1 



T,(fc) =45{G,,[0,a;(fc)]}F,(fc) 
Tt{k)^2Q{^Gy,[0,ij{k)]- 
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„/„(fc) = 43{C/^^[0,w(fc)]}F,/^(/c) 



,[0,Lo{k)]}Mt{k) 
(71) 



The last factor in each of these equations Fp{k) [Mp{k)] 
is defined as the magnitude of the total force (torque) 
carried by the TRE solution per unit amplitude. To de- 
rive Eq. (71 ) for f3 = t we have also used that the stress 
point source corresponding to a normal torque Mout at 
the origin is specified by 



-M^S'{x)S{y) 



(72) 

which can be deduced from the source's symmetries, 
and the properties: Gxy[x,y, z,uj{k)] — Gyx[y,x, z^uj^k)], 
Gyy[x,y, z,uj{k)] — Gxx[y,x, z,uj{k)]. Note that from the 
latter we also have '^{Gxx[0,i^{k)]} = 3{Gyy[0,a;(fc)]}. 
If we now eliminate k in favor of the transverse wavevec- 
tor in the support qt and define 



G(f,a;M) 



2(1 + as)qt 
F, 



G{qtr) , 



(73) 



dimensional analysis directly implies that the imaginary 
parts of the function G and of its derivatives evaluated 
at the origin are dimensionless functions of the Poisson 
ratio for the support materiaP^. This together with 

F,{k)^EbSk Mt{k)^Ck F,/h{k) ^ Ebly/xk^ 

and Eqs. (71 1, (73 1 yields scalings with kd for the trans- 
mission coefficients T^(fc) that for all branches are con- 
sistent with Eqs. ( 63 1 providing a non-trivial check of the 
reciprocity relation ( 56 1 . We have also corroborated the 
prefactor for torsion and bending in the cases where the 
aforementioned symmetries are met, by explicitly calcu- 
lating 3{G,,(0)} and 3{gf^G,,(0)-^G,,(0)}/2 and 
comparing with 



2^^ dxdy 



GyxiO) 
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dy- 



GxxiO)} 



1 

12^ 
Q{Gxx{0,a)} 



~(A) 
"t,SH 



(75) 



(59) which can be obtained from Eqs. ([60|, (|62|, ( |63| , ( |71 



( 73 1 , ( 74 1 , and the expression for qt in terms of k for 
each branch. The explicit derivation of 3{Ga;2;(0)} from 



its definition ( 73 1 is given in App. |E| We have checked nu- 
merically that the corresponding expressions for it^/ft (a) 
given, respectively, by Eqs. (E14l and (|62| coincide for all 



physical values of the sound speed's ratio a corresponding 
to positive Poisson ratios (0 < (Ts < 1/2). The analogous 
derivation for torsion (applicable to a cross section sym- 
metric under rotations by 7r/2) is greatly simplified by 
the fact that only SH waves contribute — leading to a 
universal prefactor. Finally, we note that Eq. (75 1 and 



its analog for compression allow us to interpret the u as 
dimensionless displacements and angles. 



V. RESONATOR - BATH REPRESENTATION 
AND SPECTRAL DENSITIES 

A. Resonator modes, scattering modes and their 
overlaps 

We return now to the analysis of the normal modes 
of the whole structure Uq{r) (scattering modes) and the 
resonator mode 'S'^(r) inside the beam for the purpose of 
determining their overlaps to lowest order in d/L. The 
prescription for the resonator mode given in Subsec. |IIB| 
implies that the procedures to be followed for the bridge 
and cantilever geometries differ only in the boundary con- 
ditions: clamped boundary conditions at the junctions, 
but free boundary conditions at the end of the cantilever. 
Of course the resulting Sturm-Liouville problem defines 
an infinite set of resonator modes so that ujr —^ uJn,p 
with n = 0, 1, . . . and (3 = c^t, v, h. As we are interested 
in the regime d/L <C 1, a natural requirement to identify 
these localized modes with the physical resonances of the 
whole structure is Ci!n,i3 ^n,i3 for d/L ^ 0, where the 
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former are the real parts of the poles of the propagator 
for the displacement field (18 1 — this will be borne out 



below. In this respect, it is worth noting once again that 
our problem can be viewed as a mechanical lossy Fabry- 
Pcrot (cf. Subsec. II A and Sec. Ill I. As already discussed, 

we focus on low frequencies so that ujp[k{q)], uJn,i3 ^ ^^^^ 
with ^ -k jd. In this regime the localized resonator 

modes can be associated with effective one-dimensional 
fields 4>n^p{z) defined in complete analogy to the fields 
( /)fc(q ) g(z) describing the scattering modes [cf. Eqs. (33 1, 



(35 1 and Subsec. B 1 . As expected, the refiection symme- 
tries imply that resonator modes associated with a given 
branch /3 are orthogonal to scattering modes character- 
ized by /3' ^ /3. 

We now consider the results for the reflection am- 



plitudes r(;')(0), r%,^{SS\ 



c/t 



(0), 



V /h,Sri 



(0) and for 



the ratio 6(0) obtained in the previous Section (cf. Sub- 
sec. IV A I together with the analysis for the propagation 



of low frequency modes inside the beam performed in 



Sec. Ill [cf. Eqs. (34), (36l]. It follows that to zeroth or- 



der in d/ L the scattering effective one-dimensional field 
(t>k(q),p{z) and its derivatives satisfy the same homoge- 
neous boundary conditions as those defining the (f^nsiz) 
except for the value of the field at the left junction ■ 
specified by 



0fcx/t(O) = 1 (/'fc,t,//i(0) = 1 + i . 



(76) 



On the other hand all of these effective fields are solutions 
of the TRE equation 



with 



(-1) 



(77) 



(78) 



It is then simple to use the aforementioned boundary 
conditions and Eqs. (77 1, (ITSI to obtain 



+ (-ir5^^S^^^^^I-o (79) 



where we have used integration by parts and that the 
(j)n,f3{z) are real (cf. Subsec. II B I to establish the last 
equality. We have calculated the resonator modes 0„ ^(z) 
(cf. 



Ref. |9l46j ) normalized them to the length L and 
computed their necessary derivatives [note the tt phase 
freedom for the choice of (/i„^^(z)]. These together with 
Eqs. (76 1, (79 1 and the dispersion relations for the differ- 



ent branches, yield 



,,/3|0fe(w),/3) — a/2C„,;3 ^ 



P/3 



C, t 

V, h 



(80) 

Here n = 0, 1, . . ., C^.c/t = 1, Cn,v/h = (tanh^ , 
and is the exponent of the corresponding dispersion 
relation — i.e. p^/t = 1 and Pv/h = 2. 



Equation ( 80 1 will prove useful below when using the 



overlaps to calculate the force spectral densities I{uj)n,i3- 
Naturally, its above derivation is invalid whenever k is 
a zero of the resolvent of the linear system — Eq. ( 34 1 
[Eq. (36l] for fi = c^t \(3 = v^h] — that determines the 
scattering mode (f>k,i3{z) which for k = kn (at a generic 
value of z) diverges as d/L 0. The behavior of the 
overlaps and the spectral densities in the neighborhood of 
these special points, which correspond to the resonances, 
is discussed further in App. |F] Naturally the divergent 
behavior at a;„ can be used to prove a>„ ^ w„. 



B. Relations between I{io) and p{uj), t{uj) 

In principle the overlaps calculated in the previous 
Subsection would allow to obtain the leading contribu- 
tion in d/L to the environmental spectral densities I{uj) 
from Eqs. (p6|) and (32). However the latter exact expres- 



sion has the drawback that even for a generic value of uj 
the dispersive contribution As{uj) brings into play both 
the behavior at ujr and at high frequency of the function 
S{uj), whilst the analysis we have done of the scattering 
modes Uq{f) and their overlaps with the resonator modes 
fails at these frequencies. To overcome this issue we first 



invert Eq. ( 32 1 to recover a well known expression for the 
function S{uj) in terms of /(w) 



Situ) 



u}rI{lo) 



- - tjflA/(w)] -t- iLjjiP{uj) 



(81) 



and subsequently derive from it the following approxi- 
mate relation 



{Lo + LORf {io-CoRf ^ull'^iCuR) . (82) 



Here ujr is the approximation to the renormalized fre- 
quency afforded by the solution that the equation 



Ur 



i^R. + t^R. 



(83) 



has close to ujr when \ Ai{ijJr)\/lor, <^ 1. The latter con- 
dition is necessary for the validity of the approximate 
relation ( 82 ) and the corresponding relative error is at 



most of order \IS.i{ljr)\/lor. The derivation of Eqs. (81), 
(82) is given in App. |f] where we also establish that 
\/!^i{ujr)\/u}r — > 0, I{ujr)/ujr — > for d/L and 
give the analog of the approximation ( [82| ) for the cou- 
plings {Cq}- Hence, to obtain I{uj) correct to lowest or- 
der in d/L we can replace the factor in square brackets 
in Eq. ( [82| by its limit for d/L which will be given 
by ojR — )■ ujR, I{ujr) 0. If we now consider the defini- 
tions of tq and 4'k{z), and the corresponding expressions 
for the displacement field given in Sec. |III| and substi- 
tute into the definition ( 26 1 we obtain for the remaining 
factor 



S{ui) _ S\{(l)R\(f>k(^))\'^ pf, 



LOR 



\tq\H{u^Uq), (84) 
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where we have used the definition p9|) for the overlap 



R\Uqh 



the expressions for u{f) in terms of (/'(-z) given 
in App. [B] and neglected the contributions that involve 
derivatives of (/>(^) which are higher order in d/L — note 
that (pn — > 4>B,i defined in the previous subsection, is 
normalized to the length while is normalized in the 
metric defined by Eq. (19). Here q runs only over modes 



with the branch index (3 corresponding to the appropriate 
resonator mode (frequency ujji) and we have introduced 
the number of supports 5 so that for the bridge geom- 
etry q is further restricted to right-movers. On the one 
hand the last factor on the RHS of Eq. ( 84 1 , correct to 



lowest order in kd, is just given by the appropriate P!3{uj) 
[Eq. ( 58 1 summed over all 7] which together with Eq. ( 80 ) 
for the overlaps and the relation 



Pf3 ^ » 



duj 



dk 



(w) , (85) 



valid for pp — 1,2, yields Eq. ([5|, namely 

'dujp 



dk 



2L LO 



On the other hand the last factor on the RHS of Eq. ( 84 ) 



can also be expressed exactly in terms of the transmis- 



sion coefficients Ta(w) using the reciprocity relation (56) 
derived in Subsec. IV B We can use this fact, Eq. ( 80 1 for 
the overlaps, and Eqs. (85), (84), (82), to finally obtain 
the relation (|6| between the leading contribution to I{uj) 
for each resonance and the transmission coefficients: 



/„,^(w) = ^C„,^-^(w) 



The validity of Eqs. (|5| and ([6]), to lowest order in d/L is 
affected by the caveats discussed in the previous Subsec- 
tion, when deriving the overlaps, so that the neighbor- 
hood of the other resonances is excluded and uj should be 
much smaller than top {-k /d). Note that with the judicious 
choice of dimensionless variables d/L, Lo/ujn the latter is 
not an additional approximation as Lop (tt / d) /w„ 00 for 
d/L 0. Clearly, the above procedure is questionable 
in a neighborhood of ujji where Eq. ( [80| for the overlaps 
diverges. Nonetheless, uiji is in fact included in the fre- 
quency range where Eqs. ([5| and ([6]), are valid. To resolve 
this issue it suffices to prove, without using the Eqs. (81 1, 
(82 1, that the function I{uj) remains well behaved at this 
resonant frequency as d/L (cf. App. [f|). 

In turn, it can be argued that up^n/d) always sets the 
order of the natural ultraviolet cutoff. If we consider the 
exact Eqs. (26) and (32) it is clear that the lattice con- 



stant a provides an obvious ultraviolet cutoff for the func- 
tions S{uj) and I{uj) so that the case d ^ a follows triv- 
ially. If instead d^ a, one finds that a;^(7r/d) provides a 
"soft cutoff" beyond which the functions S{uj) and I{uj) 
decay as integer power-laws. In fact for uj 3> LOp^n/d) one 
can argue that to analyze the scaling of S{uj) it is permis- 
sible to replace in its definition ( 26 1 the relevant modes 



Uq{f) inside the beam — i.e. those with Uq ^ u}0{-k /d) — 
by the normal modes the beam would present for d ^ 00. 
Thus the function ujS(u!) will scale as the DOS times the 
corresponding Fourier transform of the resonator modd^. 
The latter yields a factor of 1 /lj^ for the longitudinal di- 
rection and a factor of l/cj for each transverse direction. 
Hence we obtain S{ijj) ~ l/w™, where m = 7 for the case 
of a "bulk cross section" and m = 6 for the case of a 
SWNT modeled as a shelP^. A straightforward analy- 
sis yields in both cases 1 /uP' for the leading term of the 
asymptotic expansion of l^g{uj). The latter behavior to- 
gether with Eq. (32) leads to I{lo) ^ X/lJ with I — 3 
and 1 = 2, respectively. Finally, it is worth noting that 
the scaling for uj ^ can also be analyzed without re- 
sorting to the approximation ( 82 1 . In fact one always 



has As (0) < which together with Eqs. (|84|, (pO|, (|58l 
and (32) imply that the functions S{uj) and I(uj) share 



the same scaling as a positive power-law detailed in the 
introduction (Subsec. IA| for the different branches. 



VI. CONCLUSIONS 

In summary, we have analyzed a generic beam geom- 
etry suspended from structureless supports in the limit 
of large aspect ratio d/L ^ and provided for each of 
its low frequency resonances a Caldeira-Leggett model 
adequate to describe the associated quantum Brownian 
motion induced by phonon tunneling losses. The cor- 
responding effective Hamiltonian for the low frequency 
vibrational degrees of freedom is derived from the un- 
derlying microscopic physics performing a controlled ap- 
proximation in the natural small parameter d/L. This 
yields the lowest order contribution in the aspect ratio 
to the associated environmental force spectral densities. 
We find two general formulas for these functions Eq. ^ 
and Eq. ^ that involve, respectively, the transmission 
coefficient at a single junction and an effective environ- 
mental DOS. Whence providing two alternative pictures 
for this dissipation mechanism in terms, respectively, of 
phonon tunneling losses and support-induced modifica- 
tion of the DOS [cf. Eqs. and ([58|]. These 
yield for the spectral densities functions of frequency that 
only depend on the length L and on properties of the 
"decoupled" support and of the decoupled infinite beam. 
Furthermore they are universal in the specific sense that 
they only depend on the properties of the beam through 
the quantities that determine its low frequency effective 
theory known as TRE. 

These environmental spectral densities result in fun- 
damental limits for the intrinsic dissipation (i.e. upper 
bounds for the Q-values) with structureless supports 
which are relevant for state of the art mechanical res- 
onators in the 0.1 — IGHz frequency range (cf. Table |l] 
and Ref. [3|). It is important to note the fast degrada- 
tion as the length is shortened and that the fundamen- 
tal limit of torsional resonators greatly exceeds that of 
flexural ones (for comparable frequencies^. In fact, for 
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torsion this dissipation mechanism is likely to be neg- 
ligible when compared with other contributions to the 
mechanical damping. In addition we find that semicon- 
ducting SWNTs are far more resilient to this dissipation 
mechanism than semiconductor heterostructure realiza- 
tions (typical values for the figure of merit ksT/hujj^Q 
are at least an order of magnitude larger) . 
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In the case of 3D supports the environmental spectral 
densities are ohmic for flexural and compressional reso- 
nances and superohmic for torsional resonances, while in 
the case of 2D slab supports they yield 1// noise with 
an infrared cutoff provided by the size of the slab. It 
is worth noting that this type of noise is normally as- 
sociate d wit h interactions involving charge degrees of 
freedorrPMI! while here it arises in an insulator from 
purely vibrational effects. 

Naturally, supports can be engineered to suppress the 
phonon tunneling losses. However in many feasible al- 
ternatives analysis of phonon propagation in a beam 
geometry is essential to determine the improved limits 
attainabl^^. In this context and also when analyzing 
situations where the resonator mode couples to a pseu- 
dospin the resonator-bath representation given in this 
paper may need to be complemented with an approxi- 
mation for the complete phonon propagator inside the 
beam. In other words the effect of the "environment" 
cannot be "lumped" into the finite mechanical Q-value 
and there may be interference effects between different 
resonances arising from correlations between their effec- 
tive noise sources. Clearly, away from the resonances an 
adequate approximation for the complete phonon propa- 
gator is afforded by the lowest order contributions in d/L 
to the scattering modes [cf. Eq ( 18 )]. The latter are spec- 
ified by Eqs. (33)-(36l and our approximations for the 



reflection and transmission amplitudes given in Sec. IV 
— of particular interest is the behavior that results for 
w ^ 0. On the other hand a satisfactory solution for 
all frequencies requires going beyond the lowest order for 
the reflection amplitudes — a matter pursued in detail 
elsewherd^. 

Furthermore, we highlight the relevance of the precise 
connection, given in App. |B] between the effective one 
dimensional TRE description and the underlying trans- 
verse mode profile to scenarios where the resonator is ma- 
nipulated by coupling to an embedded optical emitteil^^l. 

Finally, we have provided a general expression for the 
spectral density of a given resonance in terms of the rel- 
evant elastic modes [cf. Eq. (10)] that allows to extend 
the treatment to other geometries. We have illustrated 
this for an axially symmetric resonator supported by a 
vertical pedestal of length at least comparable to the res- 
onator's size. A scenario which is relevant for optome- 
chanical systems based on microtoroids, microdisks or 
microspheres. 



APPENDIX A: NORMAL MODES OF THE FREE 
ELASTIC HALF-SPACE 

The free elastic half-space supports four types of nor- 
mal modes that can be classified based on the charac- 
ter of the incident wave: longitudinal bulk waves, trans- 
verse bulk waves with two polar izations (SV and SH) , 
and Rayleigh surface waveg^^ESl. In all four cases the 
corresponding eigenfunctions can be written in the gen- 
eral form: 



1 



(27r) 



(Al) 

where is the dimensionality [d^ — 2 for 7 = 5 and 
d^ — i for 7 7^ s) , Eq and q correspond to the polarization 
and wavevector of the incident wave, and A^/l, Et/i and 
qt/i to the amplitude, polarization and wave vectors of the 
refiected transverse/longitudinal waves. These depend 
on q and 7 and are given by 



longitudinal waves: 

q — p{sm cos f, sin 6 sin ip, cos 6) 
qi = p{sm 9 cos ip, sin 9 sin (p, — cos t 



qt — P \ sin 9 cos ip, sin 9 sin (yS, — y l/a — sin 



q - qi 

£0 — ^5 — ^ 1 
P P 



Al 



At = 



(A2a) 
(A2b) 



(A2c) 
(A2d) 



^asmf I , (A2e) 

2^3/2 sin 9 Vl - asiT?9 sin 26* - ( 1 - 2a siv? 9) ^ 

2a3/2sin6lv'l - asin^ 6'sin26' + (1 - 2asm^ 9)^ ' 

(A2f) 

2Y^sin 26'(2asin2 6l- 1) 

2a3/2 sin 9 ^1 - asm^9 sin 261 + ( 1 - 2a sin^ 9) 2 ' 

(A2g) 



Et = (^— COS ip\/ 1 — a sin^ 9, — sin p>\/ 1 — asin^ 



SV waves: 

q — p(sm 9 cos (p, sin 9 sin p, cos ( 



(A3a) 



qi = p [ sm 9 cos p, sin 9 sin p, — \/ a — sin^ 9 + ie 



(A3b) 

qt = p(sin 9 cos p, sin 9 sin p, — cos 9) , (A3c) 
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Eq = (cos 6 cos (f, COS 9 sin ip, — sin 0) 

= 7= > 

P\/a 

Et — — {cos 9 cos If , cos 9 simp, sm 6) , 



Ai = 



ie sin 9 sin 29 — cos^ 29 



2 \/a- sin^ 6* + sin 6* sin 26* + cos^ 20 
2\/asin26lcos26' 



2 \/a - sin^ 9 + ie sin 61 sin 261 + cos^ 26* 
SH waves: 

q — p{sm 9 cos ip, sin 9 sin ip, cos 9) , 

= p(sin cos sin 9 sin — cos 9) , 
£o = = (sin p, - cos p, 0) , 
= , Ai^Q, At^l; 

surface waves: 

q — p(cos p}, sin p) 
qi — P ^cos (y9, sin p, -~i\J\ — 

■( 

£0 = 0, e; 

£* = 7 - C^cos iv^l - C^sin v?, 1 



gt = p ( cos p, sm -iv^l - 



(A3d) 
(A3e) 
(A3f) 

(A3g) 
(A3h) 



(A4a) 
(A4b) 
(A4c) 
(A4d) 



(A5a) 
(A5b) 

(A5c) 

(A5d) 

(A5e) 



A, = 



-Cy^pCia), At ^ z^^p^/LJ=C{a) , 

(A5f) 



Cia) 



2-a«) 

[l-a«)]^/^ 
2 - ag^(Q) 



[^V)/4 + C»~l] 



-1/2 



(A5g) 



Here we have defined \/a = Ct/ci, where q and q are 
the velocities of propagation of transverse and longitudi- 
nal waves in the elastic medium, and adopted spherical 
coordinates for the wavevector q (cf. Fig. [T]). The pa- 
rameter ^ is the ratio of the velocity of propagation for 
surface waves to Ct and is a function of a that is always 
less than unit}®'. We note that for SV waves with polar 
angles 9 > cos~^ ^/l — a the longitudinal component is 
evanescent. 



APPENDIX B: NORMAL MODES OF THE BEAM 

1. TRE solutions 

At low frequencies there are four branches of propa- 
gating modes in an infinite beam: two bending branches 



(vertical and horizontal) , a torsional and a compressional 
branch. An heuristic way of understanding how these 
four types of motion arise is to decompose the associated 
displacements of the points of each cross section into an 
overall translation, an overall rotation, and a "residual" 
deformation. The vicinity of free boundary conditions 
and the fact that the wavelength 27r/fc of these low fre- 
quency modes is much longer than d determine that the 
deformation is higher order in the small parameter kd, 
when compared with the rotation and translation. Thus 
each cross section can be seen as a slightly deformed rigid 
body. Clearly if we considered a chain of coupled rigid 
bodies, each unit would have six degrees of freedom that 
would lead to six branches. However only four of them 
will have a vanishing dispersion relation and respect the 
requirement that the cross sections remain only slightly 
deformed as they move together, namely, those associ- 
ated with the three possible translations and to rota- 
tion around the beam's axes. As these considerations 
are quite general, though throughout this Section we fo- 
cus on the case when the low frequency effective theory 
can be derived from 3D elasticity [cf. Eq. Q], in the case 
of nanowires and SWNTs the expressions for the TRE 
modes that follow are warranted up to the highest order 
for which the strain vanishes — i.e. first order for bend- 
ing and zeroth order for compression and torsion — with 
discrete coordinates x, y and up to a pref actor. 

We turn now to the analysis of the small kd behav- 
ior of the corresponding normalized eigenmodes for these 
four propagating branches (/3 — c,t,v,h with m — 0). 

We may focus on a right mover given by \v^^\k)) = 
Afjfi[x,y,d,—ik]e^''^ with eigenvalue ui0Q{k,d) and ob- 
tain the left movers by reflection symmetry and (in the 
case of bending) the large decay length exponentials 
Av/h,i by analytic continuation. The problem of deter- 
mining the Taylor expansions of A/s^ and uj'pQ{k,d) in 
powers of k can be formulated as the search for an har- 
monic soluti on ug (f, t) of the 3D elastic wave equations 
for the bearrPSMl via the ansatz: 



up{r,t) 



n=0 



dz 



(Bl) 



with 4>{z,t) = e^*^^-"*). If one neglects in Eq. (|Bl| all 
orders that scale at least as the inertia (i.e. as uP)one ob- 
tains the TRE low frequency effective theorjEllll. Then 
the function (t>{z,t) will correspond to the effective one 
dimensional field. On the other hand, the lowest order 
for Lu"^ yields the corresponding approximate dispersion 
relation lo cx fc^'' . It is straightforward to realize by sub- 



stitution ofEq. (Bl) into the elastic wave equations for 
the bearrP^I^ that the result of this truncation at order 
2pp — 1 will yield an exact solution of the corresponding 

static equations (w 0) provided that vanishes. 
Thus the leading terms of the Taylor expansions of the 
transverse profiles Ap Q can be extracted from the corre- 
sponding static solutions. These are given by 
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compression: 



Ux{x,y,z) 



-aa:-(z), 



Uy{x,y,z) 
Uz{x,y,z) = (p{z) 



oz 



(B2) 



torsion: 



Ux{x,y,z) 
Uy{x,y,z) 



-y(t>{z) , 
x(t){z) , 



Uz{x,y,z) = tl^{x,y) — {z) ; 



(B3) 



vertical bending: 



Ux{x,y,z) = 0(z) + -(a;^ - y^)^(z) , 
2 az^ 

Uy{x,y,z} = axy—^{z) , 
oz^ 

d(f> d'^cf) 
Uz{x,y,z) = -x — {z) + axv{x,y)-^{z) . 



(B4) 



In the above, the coordinate system is oriented along the 
principal axes of inertia, {x,y) = corresponds to the 
center of mass of the cross section, and a is the Poisson 
ratio of the beam's material. The functions ip{x, y) and 
Xv{x, y) are determined by the static equation V • <t = 0, 
the free boundary conditions, and the condition {uz)s — 
0. In the case of bending these yield 



dXv, ^ dxy xy 
— ia:,y)=—{x,y) + -, 

{Xv{x,y))s = 0, 



y 



2a 



(B5) 



where Xv is an harmonic function that at the cross sec- 
tion's boundary (a;[^],y[Z]) satisfies 



1 



X.(0=^[3(l + a).T2y + (l-a)y^]„ 



2(1 + 0-) 



dx 



y—dl'. 



(B6) 



The analogous relations for the case of torsion are given 
in Reference '45] while the transverse profile of the dis- 
placement field for horizontal bending can be obtained 
from the RHS of Eq. (|B4| and Eqs. p5]), (|B6]) via the 
ansatz x ^ y, Xv{y,x) -> Xh{x,y), Xvi^x) Xh{x,y). 

Finally, the desired approximation for A0 q[x, y, d, —ik] 
is obtained from Eqs. (|B2|-(B4| via the replacement 



dz" 



■;^=^(zfc)" for (3 = c,v,h 
j^^iikr ior(3 = t 



(B7) 



where we have taken into account that the ID continuum 
modes \v''^\k)) are normalized in the standard Euclidean 
metric. Whence, we define 



/3 (2) 



|^M,,,(0,0,z) for/3 = c 

^ z-V xug{0,0,z) for/3 = i 



\/2t, 

^Ux/y,^iO,0,z) 



for f3 — v/h. 

(B8) 

The corresponding approximate dispersion relations Q 
are specified by 



/ C 
Pbh 



' Ebly/x 

PbS 



(B9) 



Here li are the moments of inertia with respect to the 
principal axes of a cross section per unit surface density, 
C is the torsional rigidity of the beam and Eij its Young's 
modulus. 



2. End corrections 

If one considers now a semi-infinite beam (z > 0) it is 
clear that complex values for k (the wavevector along z) 
are physically meaningful provided that their imaginary 
part is positive. The corresponding solutions can be un- 
derstood as the analytic continuation of the "traveling 
wave" eigenmodes of an indefinite beam. More precisely, 
Lu'^{k,d) and the transverse profile A{x,y,k,d) will be 
multivalued analytic functions of k that yield harmonic 
solutions on the paths in the complex plane specified by 
3?[tj2(/c,d)] > 0,3[w2(fc,d)] = 0. This can be viewed 
explicitly for a cylindrical cross section given that the 
resulting problem is separable (see Ref. [Hj and refer- 
ences therein) but will hold in general provided that the 
boundary of the cross section is well behaved. We as- 
sume for simplicity that the latter is characterized by a 
single length scale d but all our considerations follow for 
an arbitrary cross section (within the caveats that fol- 
low) provided that we reinterpret d as the largest chord. 
Heuristic considerations imply that for a generic geome- 
try the finite complex zeroes of (^^(fc, d) will have a real 
part that is at least of order 1 /d. We will further assume 
that at these zeroes w^(fc,d) is non-singular (i.e. they are 
not branch points) and has non-zero derivative; and that 
the associated path yielding physical solutions admits a 
parametrization in terms of ^[k] or 3[fc] that can be Tay- 
lor expanded. Thus each of these non-trivial zeros yields 
a branch of end corrections (with label^l^ /?, m as de- 
fined in Subsec. IV A I specified by the transverse profiles 
Ap^mix, y, d, Kp^mi^'', d)] and the "wavevectors along z" 



Hf3^m{^ T d) = ~ik, where the functions j4^.,„, K^.m can 
be expanded in natural powers of w^. 

To analyze the problem of transmission at a 3D- ID 
junction it proves useful to consider the overall displace- 
ment and rotation of a given cross section. The former is 
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simply defined as the displacement of the cross section's 
center of mass. In turn, a natural definition of a spa- 
tially averaged angle for a neighborhood V undergoing 
harmonic motion is afforded by considering the maxi- 
mum total angular momentum Ly with respect to the 
center of mass fy of V over a period: 



1 



(u(f))ang,V = - I '^Lv 
UJ 



(BIO) 



where I is the equilibrium inertia tensor for V . If we 
orient the axes along the principal axes of inertia and set 
the origin at fy, we obtain: 



{U{f))s 



dr^ r X u{r) ■ i 
/ydr3(r2-xf) 



(Bll) 



If V is small enough this reduces to 

(w(^))ang,y • I w X u(0) • i + jiUjk{0) with eikj = 1 

(B12) 

where 7^ = eijk{Ik ~ the relative order 

of the corrections is at most u{Q)d/ min{|uij(0)|} with d 
the typical dimension of the neighborhood V . We note 
that throughout Sec. llVlwe use the notation 



~e^{u{f)) 



ang,S 



(B13) 



where S is the beam's cross section at the origin. 

In the limit w — s- each branch of end corrections yields 
an exponentially decaying solution (cx e"'''' '"^^''*'^) of the 
corresponding static problem (i.e. a static end correc- 
tion). To establish Eqs. (46 1 underpinning our analysis 
of a 3D-1D junction, in addition to the aforementioned 
analyticity properties, we just need to use the following 
exact universal properties of static end corrections: 



.771 

(0,d)]e 



}>5= 

(B14) 

(F • [Afi^^ [x, y, d, «:^,™(0, d)] e-'='^-(°''')^})a„g,s = 

(B15) 

(A^,„i [x, y, d, K/3,„(0, d)])s = 

(B16) 

(A/3.™ [a;,y,d, K/3,„(0,d)])ang,5 = 

(B17) 

{[ dr'' G (r - f', 0) • F- {A/3,„[a;',y',d,At^,,„(0,d)] 



X e 



-K^,„(0,(i)z' _ Q 

(B18) 



(/ dr'' G (r - f', 0) • Y-{Ap,m[x\y\d,Kp^m{Q.d)] 
Js 



X e 



})ang,S = , 
(B19) 



where now G(r, 0) is the static Green's function of a free 
elastic half-space loaded at the free surface. Equations 



(B14l and (B15I correspond, respectively, to the total 



force and torque applied at S and can be established 
by simply stating the equilibrium conditions for a given 
finite segment of the beam. In fact the exponential de- 
pendence on z implies that the total external force and 
total external torque applied to the finite segment result, 
respectively, proportional to the L.H.S of Eqs. (B14) and 



(B15I 



On the other hand Equations (BI6I and (B17l can be 
reduced to the following Lemma: 

Lemma: For a finite beam of length L in equilibrium 
subject to specified displacements at z = and free bound- 
ary conditions at z = L the average displacements and 
angles at both ends coincide, 

by taking the displacements at 2: = specified by 
A^^m[a;, y, d, K^^m(0, rf)] and then sending L — > 00. In 
turn, the above Lemma can be understood using a vari- 
ational argument given that the elastic energy density 
is positive definite in the relative distances between the 
material points of the beam (see Ref . a nd references 
therein). Equations (B14l through (B17) underly the 



standard recipes for the boundary conditions in the TRE 
treatment of static small defiections (i.e. "linearized- 
strain" theory j^^I^. In this respect, it is important to 
note that for a beam in equilibrium the corresponding 
TRE solutions, already discussed, are exact within the 
linearized-strain 3D theory. Thus the validity of the lat- 
ter recipes hinges on whether the strain can be linearized 
and not on the smallness oi d/L. 

We now consider the cantilever geometry (cf. Fig. [ij 
for an arbitrary length L and apply the stress source as- 
sociated to a given static end correction (/3, m) dX z — L, 
the usual free boundary conditions at the other surfaces 
and fixed displacement boundary conditions for z — > —00 
in the support (the latter ensues for w — s- given the 
boundary conditions satisfied by the support's Green's 
function). Equations (|B14|)-(|B15[) imply that the TRE 



part of the corresponding static solution inside the beam 
[cf. Eq. (B2|-(B4)] is at most an overall displacement and 
rotation. However a variational argument analogous to 
the one underlying the above Lemma implies that this 
displacement and rotation of the beam should vanish so 
that the solution is given solely by "end corrections" . 
Then, continuity at the junction yields Eqs. (B18l-(B19l 
if we take the limit L ^ 0. We note that this derivation 
will be valid for any well behaved support irrespective of 
any symmetries. 



APPENDIX C: GENERALIZATIONS TO 
ASYMMETRIC STRUCTURES, 2D SLAB 
SUPPORTS, AND SWNTS 



The derivations in Subsec. [IV A| rely on the assumption 
that both the support and the beam are symmetric under 
refiections with respect to the x and y axes. It is straight- 
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forward to extend these results to the case where only the 
beam's cross section presents these symmetries so long 
as we assume that the (now asymmetric) support is still 
characterized by a structureless continuum of free eigen- 
modes Mq""* (f) that has an effective 3D density of state^^. 
More specifically each mode (f) should be character- 
ized by a single length scale of order 1/qt (where qt{oJq) 
is the wavevector for propagation of transverse elastic 
waves in the support's material). Thus the analogous 
property will hold for the corresponding Green's func- 
tion G(f , f', w) (f' lies at the boundary where the normal with 



stress vanishes when the support is unloaded). A simple 
example would be an elastic quarter-space or any fraction 
of an elastic space with origin at the junction subtending 
a solid angle of order 7r/2. In this more general scenario 
there will be mode mixing between the branches due to 
the scattering at the support. Equation (38 1 will still 



be valid with the analogous definition for but now in 
Eq. ( 39 ) the decompositions of the fields Am+ and in 



terms of the eigenmodes of the beam will involve an ad- 
ditional sum over all the branches (3 (i.e. — *■ ^). 
Concomitantly these fields will no longer have a branch 
index except for in case (i) where (3' will denote the 
character of the associated incident field Win,^', and the 



analog of Eq. ( 45 ) will no longer involve a symmetriza- 
tion operator Sp. The support mode labels q no longer 
bear any relation with the beam branch index /3 and now 
a convenient parametrization of the complete solution is 
afforded by qt instead of k. 

In complete analogy to the symmetric case we substi- 
tute the modified decompositions ( |39[ ) into the modified 
Eq. (45 1 and Eq. (38 1 and take on both sides the spatial 
averages (. . .)s andT?. .)ang,s- Note that the beam's sym- 
metries imply that on the LHS of each of the resulting 
equations for the amplitudes and cpm only the term 
with the relevant (3 — i.e. the one associated with the 
TRE branch for which the corresponding average does 
not vanish — will survive. Once again we solve for the 
TRE starred amplitudes in the linear system arising from 



the analog of Eq. ( 45 1 and substitute them into the lin- 
ear system arising from the analog of Eq. ( 38 1 , whence 



amplitude (whose scaling as qtd is preserved) arising 
from Ct,Q{qtd) can be readily added to yield 



,(0) _ '-"^t,Q 
" dq, 



J2 



(0) qt + 

{f3^t,m}(ETRE 

dr'^G(r-r-',0).z4°^(r-'))ang,s41(0) 



(CI) 



/1;L(^i>^2) = i-ir^xy,k'^^,{q,) m = 0,1 



V2^ 



/io(a;i,a;2) = i^^k^iqt) , 

V ^TTO 



En 



(C2) 



where we have used the transverse profiles of the TRE 
modes given in App. |B]and G (f, r',0) is now the static 
Green's function for the asymmetric support under con- 
sideration. 

The aforementioned results for the starred amplitudes 
imply that for incidence from the beam to lowest order 
in qtd there is no mixing between the branches so that 
the refiection amplitudes for [3 ^ j3' vanish and those for 
(3 = P' coincide with the r^^^l[k{qt)d\, r^^^\^g^[k{qt)d\ al- 
ready found for the symmetric case. On the other hand 
for incidence from the support the solution is a superpo- 
sition of all four b ranc hes with approximate amplitudes 
given by Eq. (Gil and for [3 ^ t the same results 

as in Subsec. 



IV A 



we solve for the un-starred TRE amplitudes. As before 
the end corrections yield terms that scale at most as the 
inertia (qtd)'^. Thus, we find that — with the noticeable 
exception of the un-starred torsional TRE amplitude for 

case {ii) — the results for the amplitudes ci*^^ and cp„i in the prefactor t^^L Thus, in general the results for the 



in terms of the modified Uq (f). If 
we now revise the treatment of the modes of the whole 
structure given in Sec. |III| the fact that the branches do 
not mix implies that to lowest order the displacement 
field inside the beam (and away from the junctions) is 
just given by adding coherently the four contributions of 
the TRE branches which are given by similar expressions 
as in the symmetric case but with modified transmission 
amplitudes (i.e. the reflection amplitudes and the 
ratio b remain the same). Furthermore, the symmetries 
of the beam imply that a resonator mode of branch f3 
only has non-vanishing overlap with the contribution to 
a scattering mode of the same branch so that the quan- 
tities {u'j^\uq) (to lowest order in qtd) coincide with those 
for the symmetric supports except for the modification 

(0) 



are given by Eqs. (44) and (47 1- (50) with the following 
straightforward modifications: kd is replaced by ^/qtd, 
u'q^ (f ) denotes now the free modes of the "generic" sup- 
port, and in case {i) the RHS of the equations for the 
starred TRE amplitudes are now multiplied by Sp^/ji . 

The caveat for incidence from the support [case (m)] 
is that now Ctfl{qtd) has a contribution comparable 
to c|*g'((7t(i) of order qtd arising from the other TRE 



force spectral densities (cf. Sec. VB ) that we will extract 



branches. Thus in this case Eq. (41) is no longer valid. (cf. Sec. VB and App 



However the correction to the approximate transmission 



from these overlaps will be completely analogous to those 
for the fully symmetric case. 

The only caveat is when there are degeneracies between 
resonances of different branches, i.e. Wm,/? = tJn./S' = 
with [3 ^ (3' . Unlike the fully symmetric case there is 
now mode mixing which, albeit higher order in d/L, 
may nonetheless invalidate our approximations for the 
functions Im.piuj), In.p' (t^) in a neighborhood of LUfj 
pp. pl). The situation is completely 
equivalent to the failure of our treatment for In,[3{^m.p) 
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with n m with the difference that now the syndrome 
occurs precisely where the function is more relevant. On 
the other hand it is simple to realize that this breakdown 
is completely consistent with the nature of our approxi- 
mation which is perturbativc in d/L. If the asymmetry 
removes the selection rules that prevented the supports 
from coupling the two degenerate modes, the Caldeira- 
Leggett Hamiltonian provided by Eq. 
discrete resonator mode no longer 
representation as d/L ^ 0. In particular, this will affect 
the two bending branches for cross sections symmetric 
under rotations by 7r/ 2 wh enever the supports break the 



( 20 1 with a single 
affords a convenient 



Ry symmetry [cf. Eq. (14 1 and Fig. 111. Thus, in a realis 



fmi 



'-y 

tic scenario (i.e. at finite d/L and finite mode-splitting) 

our results for the bending resonances 0Jn,v/h will only 
be applicable to the case of SWNTs and nanowires (dis- 
cussed below) provided the mode-splitting is negligible 
compared with the natural linewidth oJn.v/h/Qn.v/h in- 
duced by the phonon tunneling. 

We have also analyzed the case when the cross sec- 
tion is also asymmetric following the analogous proce- 
dure. The key point in this further extension is that 
the lowest order contributions in kd to the TRE mode 
profiles given in App. [B] satisfy the reflection symmetries 
for arbitrary cross section (with the axes oriented along 
the principal axes of inertia). Thus when the starred 
amplitudes are zeroth order in qtd the overlaps 
will once again coincide [to lowest order in qtd and up 
to a modification in the prefactor with those for 
the symmetric supports. It follows that the extension to 
asymmetric cross sections carries over for the results that 
will ensue for the force spectral densities (cf. Sec. VB) 
except for torsional resonances — in which case the rele- 
vant transmission amplitudes are higher order — and the 
aforementioned case of degeneracies. 

Furthermore, we have used the same framework to an- 
alyze the analogous problem of an abrupt junction be- 
tween a rectangular beam of width w and a slab (sup- 
port) of the same thicknesJ^Hl t ^ w = c?. Thus, we can 
use instead of 3D elasticity th e 2D effective theory pro- 
vided by "thin plate elasticity"'22llll adequate for phonon 
wavelengths much larger than t. A completely analogous 
treatment is feasible based on the 2D analog of Eqs. (38 1, 
(39 1, with: (i) ulj'\y,z) and G(?/,z,a;) denoting now, re- 



spectively, the free modes and retarded Green's function 
of the two dimensional elastic half-plane (i.e. a semi- 
infinite thin plate lying on the j/z-plane in our notation), 
and suitable redefinitions of {ii) the operator F and {Hi) 
the beam's mode profiles Afj „i[y, d, Kp^m{k, d)] which are 
now one-dimensional. The reflection symmetries also ap- 
ply in this case so that the vertical bending and torsional 
branches of the beam only couple to the flexural modes of 
the plate while its horizontal bending and compressional 
branches only couple to the in-plane modes of the latter. 
Thus one can prove the same results for the reflection 
amplitudes and obtain results analogous to Eqs. (47 1, 
(|50| for the transmission amplitudes tq, htq in terms of 
the averaged displacements and angles of the slab sup- 



port — the redefinitions {iii) imply that in the corre- 
sponding cross-section-dependent prefactors one has to 
replace S w and 1^ \N{r^)s (cf. Subsec. I A). The 
only caveat is that for vertical bending the second rela- 
tion for the starred amplitudes in Eq. (|49|) is no longer 



homogeneous since now the averaged angles 9y scale as 



y/uj for the relevant "free modes" of both the beam and 
the plate (flexural modes). This leads to a modification 
of the ratio b for vertical bending that now may acquire 
a g-dependence. As a consequence of the above anal- 
ogy the relations ^ [(|6|] between the spectral densities 
In,j3{^) and the effective environmental DOS of the sup- 
port P/3(tf) [transmission coefficients t^(w)] will be the 
same (in terms of the analogous functions for the thin 
plate "support" ) as in the 3D case with the exception of 
the expression for Cn,v that will be modified for small n 
— for large n the large decay length exponentials can be 
neglected and the dependence on the ratio b of In^piuj) 
becomes negligible. 

Finally, we turn to the case of SWNTs and nanowires. 
One should note that as there is no valid underlying "mi- 
croscopic" theory applicable to both the supports and the 



beam the precise analog of Eq. ( 38 1 will in general be un 



known, depending on details of the clamping procedure. 
However, one can argue that insofar as there is no dissipa- 
tion inside the "junction" the above results for the reflec- 
tion and transmission amplitudes will still hold — with 
the straightforward replacements in Eqs. (47 1, (49), (50 1 



that result from the redefinition of the transverse mode 
profiles: S — > Nc, Iz Nc{r'^)s, with Nc the number 
of atoms in the unit celP^ (of the appropriate equivalent 
ID chain in the case of a nanowire^. In particular, one 
should note that the analog of the universal properties of 
the end corrections (cf. App. |b]) are expected to hold if 
the SWNT or nanowire is to be regarded as "clamped" . 



APPENDIX D: EXTENSION TO 
MICROTOROIDS AND OTHER GEOMETRIES 



The facts that the resonator mode and the scat- 
tering modes \uq) are solutions of the time-independent 
elastic wave equations (cf. Sec.|lT]), and that u'j^{r) is real, 
directly imply^Sl 



^R{u'jj}Uq) ^ 



1 



dr^ U 



da, 



Ps JVr 



1 



R,i' 



1,^3 



Ps JVr 



— / dr^Ua 



da 



R.ij 



(Dl) 
(D2) 



where we use Einstein's sum convention and that u'f({f) 
only has support in Vr. The above can be re-expressed 
as 



^l{u'ji\Uq) 



1 

Ps 



_d_ 



''R.i'^qM ) 



^R,ij"q,ij 

(D3) 



26 



1 



Ps JVi 



_d_ 
dxi 



q,i^ R,ij J 



R,ij 

(D4) 



by using that the stress tensor is symmetric and that the 
contraction of a symmetric tensor with an antisymmetric 
one vanishes. As the mapping of the strain Uij onto the 
stress (Tij corresponds to a symmetric matrix (given that 
its entries are the second derivatives of the elastic energy 
with respect to the strain components) we have 



u 



R.ij^l-ij 



u„ 



q,ij'^R,ij 



(D5) 



Substracting Eq. (D4) from Eq. (D3l and using Eq. (D5l 



the divergence theorem, and the free boundary conditions 
satisfied by the displacement fields at the free surfaces of 
the resonator we obtain 



(D6) 



which together with Eqs. (26 1 and (82 1 finally yields 
Eq. ( 10 ), namely 



s 



(u'r ■ cr, 
c^(q)] . 



^r) -n 



very general as it only assumes 



The above 

\/!^i{lor)\/ujr <C 1 to ensure the vahdity of Eq. ([82l 



We now focus on an axially symmetric structure con- 
sisting of a "resonator volume" supported by a vertical 
pedestal. Concrete relevant realizations of this geome- 
try are microtoroids (cf. Fig. [2|, microdisks and micro- 
spheres. We consider the regime in which the contact 
area S between the resonator volume and the pedestal 
satisfies S <^ hD and \/S < h, where D and h are, re- 
spectively, the largest diameter and smallest characteris- 
tic dimension of the former. Thus there will be axially 
symmetric^ resonances with typical wavevectoiJ^ ku that 
satisfy VSkfj <C 1. We focus on one of therrP^ and as- 
sume that it is "isolated" (cf. Subsec. |I A^ . The limit 
5 — > yields a finite resonator volume subject to free 
boundary conditions on its whole surface, which specifies 
the natural choice of boundary conditions for the cor- 
responding resonator mod^^ u'j^. This directly implies 
that the terms involving cr^ vanish in Eqs. (D6), (10) 
while property VSkff, <C 1 allows us to factor u'j^ out of 
the integral which reduces then to the total force across 
S associated with the scattering mode q. 

We assume that for studying the propagation of modes 
with wavevector q ^ kji the pedestal can be modeled 
as an infinite beam with slowly varying cross section 
qS{z)/S'{z) ^ 1. This should capture the adiabatic limit 
of perfect impedance-match with the substrate for the 
purpose of studying the dissipation. In turn, our prior 



treatment of a 3D-1D junction (cf. Sec. IV) can be ex- 



tended to our present context (with the roles of support 



and resonator interchanged). Now the 3D object (res- 
onator volume) presents finite dimensions comparable to 
the relevant wavevectors. Naturally, providing rigorous 
derivations (as the ones given in Sec. IV I is now compli- 
cated by the lack of an explicit expression for the Green's 
function of the 3D object. Nonetheless, heuristic con- 
siderations imply that: (i) in the limit VSq <C 1 with 
\/~S < h the scattering mode Uq can be asymptotically 
approximated in the pedestal by (which corresponds 
to an incoming wave satisfying clamped boundary con- 
ditions at S) provided Lij{q) is not close to a resonance 
LUji of the resonator volume, and (ii) once again I{uj) is 
smooth in a neighborhood of tun. The latter implies (as 
for the beam geometry) that the restriction to modes Uq 
with uj{q) ^ LOn is immaterial as it only invalidates our 
result for /(w) at lu = lj^ ^R- In turn, it is straightfor- 
ward to extend the formalism in App. [B]to a beam with 
adiabatically varying cross section. Thus the universal 
properties of the end corrections imply that to lowest or- 
der in %/ Sq only the TRE part of Uq contributes to the 
aforementioned total force. The axial symmetry of the 
resonance further implies that only compressional modes 
yield a finite contribution^^. Thus we obtain (cf. Fig.[2| 



dr^ u' 



R 



-E.\l^_u'i^{0)-z^{0) (D7) 



,d(j)q 

dz 



where the corresponding effective one dimensional field 
satisfie^Sl 
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^'^'^' + f-^(9)<^, = (D8) 



S dz 



and the approximation by m*^^ implies <j)q{Q) = 0, 

— 1 — here we define as the in- 

coming (outgoing) components of the scattering mod^^ 
q. We adopt for simplicity an exponential dependence 
S{z) = 5*6"^^ which leads to the following 



Ct>q{z) = 4>[-\z) ~ 4>'^ + \z) . 



(D9) 



with 



dq 



£2 

4 



(DIO) 



Then, substitution of Eqs. (D7l, (D9l, (DIO I into 



Eq. ( 10 ) allows us to obtain 




ujjQ{uj~uJi) . (DU) 



Clearly the infrared cutoff uji — \J E^j pj^ jl is an arti- 
fact of the adiabatic assumption qS{z) / S' {z) ^ 1 which 
implies a; 3> w/. In the latter appropriate limit we 
finally obtain Eq. for the spectral density where 

wr,z(0) = ^J mji/ psu'pj 2(0) corresponds to the resonator 
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mode normalized to the "relevant coordinate"!^ (in our 
specific context AZ5/2). Note that, as expected, the re- 
sult is independent of F. As for the beam geometry, for 
typical materials the Q will be mostly a size-independent 
geometric property. 



APPENDIX E: DERIVATION OF 9{G^:j,(0)} FOR 
THE ELASTIC HALF-SPACE 

In the following we define u = G ■ x and use dimen- 
sionless variables setting qt = l. We need to solve for the 
half space {z < 0): 



-V(V-u) + u = 0, (El) 



1 -2a, 

subject to the following boundary conditions at z = 0: 

2u^, = Six)Siy) , (E2) 
2 

n [(1 ~ <^s)Uzz + <Js{Uxx + Uyy)] = . 

(E3) 

One can can construct the solution as a superposition of 
longitudinal (j = 1), SH (j — 2) and SV (j — 3) waves: 



u{x,y,z) 



47r2 



CO J —CO 



where we have defined 

4 _ I- 2a, 



(E4) 



c\ 2(1 -a,) 



a(o-s), a2 = a3 = l, (E5) 



^1 = {qx,qy,qz.i), 



£3 = {-qzMx, -qzMv^Qx + qy) ■ 



£2 = {qy, -qx,0) 



(E6) 
(E7) 



We note that q^j = — ^ aj — q'^ — q"^ can be imaginary, 

we do not need to normalize the Sj for our purposes and 
a < 1. 

In order to obtain the outgoing solution one can intro- 
duce damping, then calculate (E4| that will correspond 



to the steady state solution, and, finally, take the limit of 
damping coefficient going to zero. If we decompose the 
total displacement field into transverse and longitudinal 
components u = u± + U|| , the modified equations of 3D 
elasticity can be written as 



+ e — - 

dt^ dt 

V X MM = , 



■q 



V • uj_ = , 



(E8) 
(E9) 



with e > 0, and cj^ = qtCt, c\\ — qtci. We look for solu- 
tions u{f,t) = eje*('j''''^"*\ which leads to: 



1, 



1 



e 

i — 



(ElO) 



with ci = C|| and 02,3 = c^. These solutions can be 
obtained from the undamped ones by the replacement 
ct.j aj(l + it). The analysis of the limiting procedure 
e ^ allows the determination of the adequate integra- 
tion contour C in the complex plane in the standard way. 



If we substitute expression ( E4 1 and Fourier transform 
we arrive at: 



qx{ql + qy- ql.2) 



2q 

2%gza -qxqz,2 qy{ql + ql- qla) 

{l-as)qli + a,Ml + ql) {\ - 2as)qza{ql + qD , 



We adopt polar variables qx = p cos if, qy — psimp and solve for the Uj 

^ P\J\ — p"^ + ie cos Lp 



ui{p,ip) = i{l - 2c7s)- 



D{p) 



sm ip 



where 



Dip) =2(1 - 2as) 



. [(1 - as)a{cTs) + {2as - l)p^] cos<^ 



pD{p) 



(Ell) 



(E12) 



(E13) 



If we consider the integral of expression (E4| evaluated at z = over the contour C we can identify a regular and a 
singular contribution to the displacement at the free boundary. We are only interested in the imaginary part of the 



28 



X component to which only the regular part contributes. Finally, evaluation at the origin x = y = Q leads to 



TT 



where we have eliminated is in favor of a, the contri- 
bution of the Rayleigh pole associated with the surface 
waves (SAW) is given by 



1 {i-ei'i?Vi^+{i-e)VT^ 



16^ 6(a-l) + 2(3-2a)$2_^4 



(E15) 



and we have defined 

p{v) = 16(a - l)v^ + 8(3 - 2a)t;2 - 81; + 1 . (E16) 

We note that the parameter ^, already introduced in 
App. [a] satisfies ^ = where is the only real 

root of ( |E16[ ) greater than unity. 



APPENDIX F: DERIVATION OF EQUATION (82| 



We will first invert Eq. (32 1 to obtain Eq. (81). To this 



effect we define a complex function G{uj) on the real axis 
such that for > 

3[G(w)] = -5(w) and 5R[G(w)] = As(w) (Fl) 



and the extension to a; < is specified by 



(F2) 



which, after rearranging the integral for l^s{^) using the 
partial fraction expansion of 2uj' /{u}"^ — co''^), yields for 
the whole real axis 



dcj I sgn(cj)i)(|w|). 



Given that S{lo) is well behaved (Coo for lo > 0) and 
given its behavior for w ^ and — > 00 discussed at 



the end of Subsec. V B Eq. (F3 1 implies that the real and 
imaginary parts of G{u!) are Hilbert transforms of each 
other. Then, it follows from Titchmarsh's theorem^ ^ that 
G{z) is analytic in the complex upper half-plane and that 
the integral over a semicircular contour in the latter tends 
to zero as the radius is increased. Thus we arrive at 



Giz) 



1 
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UJ - z 



(F4) 



valid for 3(z) > 0. We note that G{uj) gives the propaga- 
tor of the resonator mode's canonical coordinate at zero 
temperature so that in physical terms the above prop- 
erties of G{z) can be understood as a consequence of 
causality — the inverse Fourier transform of G{uj) only 



has support for t > 0. If we now substitute Eq. (F3 1 into 



Eq. ( F4 1 , interchange the order of the integrations over 



UJ and Lu' , perform the former, and rearrange the latter 
we obtain 



G(z) = 



1 



duj'- 



2u' 



72 ' 



(F5) 



which with the help of Eq. ( 27 ) can be re-expressed as 



1 



1 



duj- 



To proceed we define the function 
E(z) = — tj^ — 



1 



S{lo) . (F6) 



(F7) 



which using Eqs. (F3l, (Fll and (32 1 can be shown to 
have the property 



-sgn(u;)^(|c.|) 
Al{u;) + SH\u;\) 



-ujR sgn(tj)/(|w|) 



(F8) 

As already discussed in Subsec. VB S{lo) has a natural 
ultraviolet cutoff, given by oj* ~ LOpij: /d), beyond which 
it is bounded by a power law l/w™ with m > 6. From 
this behavior one can derive that the last term in the 
RHS of Eq. (|F6| is bounded by \a.\z/oj^\/\z\^ if m = 6 



and by l/|z|°if to > 6. This property together with 
Eqs. (F7l and (F6I directly implies that 



\nz)\<o 



In |z/a;*| 



or |S(z)|<0 



(F9) 



(F3) respectively, for z 00. On the other hand Eq. (F5) im 



plies that the propagator G{z) has no zeros with 3(2) > 
so that S(z) is analytic in the upper half-plane. This 
property and Eq. (F9l allow us to obtain 



1 f°° y 



, S(c.') 



(FIO) 



using the Cauchy integral formula. If we now take the 
real part of this equation and use Eqs. (F7l, (F3l, (Fll 



and (F8I we obtain after some simple rearrangements 

As(w) 



A|(c.) + 52(^) 



(Fll) 



for LO > Q. If we assume that the functions 1^i{lo) and 
I{uj) are known, Eq. (Fll I together with Eq. (32l provide 



a system of algebraic equations for the unknowns l^giuj) 
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and S{uj). Thus we can solve for S{uj) to finally obtain 
the desired result given by Eg. ([sT] ). 



If we now consider Eqs. (Fll ) and (F8l, we have 



E(cij) = LUj^ [A/(cj) — i sgn(Lu)I{\L 



(F12) 



that together with Eqs. (|20[) and (21 1 allow us to inter- j^q-^ 



pret the function T,{uj)/2u;r as the resonator's self-energy 
induced by its interaction with the environment. Given 
the properties of /(cj) discussed below wc may assume 
that 5^(2:) is analytic in a neighborhood of z = LUfj. Thus, 
it follows from Eq. ( |F7| that any pole of G{z) in that 
neighborhood is given by a root zji of the equation 



9"E 



n—l 



- (l^r) {z - Luj,)" = . (F13) 



In turn, the behavior of the function S^uj) discussed be- 
low [cf. Eq. ( |F25| )] and Eq. (|F3| imply that for d/L 
the propagator G{z) has indeed a pole in the lower half- 
plane close to ujn that tends to lur. If we now consider 



the asymptotic expansion of {zr — uifi) /luh and E(aj)/a;^ 
as a function of d/L for d/L — > and substitute it into 
Eq. ( F13l, it is straightforward to realize [cf. Eqs. ( |F8[ ), 
(Fll I J that {zr — u!r)/u!r — > directly implies the prop- 
erti cs \A i{ujr)\/u!r — + 0, I{ujr)/ujr — + used in Sub- 
sec. 



V B and in the derivation of Eq. ( 82 ) pursued be 



Heuristically, this vanishing of the self-energy can 
be viewed as an unavoidable consequence of the behavior 
of the transmission coefficients r(aj) for d/L ^ un- 
veiled in Sec. II VI In fact the interaction with the bath 
giving rise to the self-energy can be understood in terms 
of phonon tunneling between the beam and its supports, 
i.e. a mechanism that is suppressed as t(u;) — *■ 0. 

To proceed we reformulate Eq. ( 82 1 by defining (uj) 
as its RHS. Thus, it is enough to prove 



I{U!) — L{'jj) 



1 + O 



(F14) 



This can also be written as the condition < 
\Aj{lur)\/u!r where the function E(lu) is defined by 



E{uj) 



Ai{uj) - Ai{CdR) - i [I{uj) - I{CdR)] 



UJ — UJfl 



Ai{Qr) 



LO — UJR + i 



(F15) 



[cf. Eq. (83)] so that given Eq. (81) we have 

2 



1 



LUr, + LU 



E{uj) 



(F16) 



We note that Eq. (83) is just the real part of Eq. (F13) 
restricted to the real axis. If we use again its asymptotic 
expansion in terms of the aspect ratio d/L it is straight- 
forward to derive from Cjr ujr when d/L ^ that 
2{u}R — ujr)/Ai{u!r) 1, which together with Eq. (83) 
implies that Aj{lur) and Aj{ujr) coincide to lowest or- 
der in the aspect ratio. It is worth noting that though 
the real part of the exact pole of G{z) differs from ojr, 
as the imaginary part of Eq. (F13 1 comes into play, anal- 
ogous considerations yield 2{zr — lor) / Ai{ljr) — s- 1 so 
that Ai{ijjR)/2 provides the lowest order contribution to 
the renormalization of the bare frequency. 

In order to ensure the validity of Eq. ( [82| ) wc find that 
in addition to \Ai{lor)\/ujr ^ 1 the following assump- 
tions (whose validity is discussed below) are needed: 



|A/(c.) 



Qnj 



UJ" 



(F17) 



for n — 1,2 and frequencies smaller than the ultravio- 
let cutoff 07* 3> LOR. Heuristic considerations imply that 
from Eq. ( |F17 1 and the behaviors of Ijuj ) for cj and 

VBI), it follows 



UJ 00 (discussed at the end of Subsec. 
that the only relevant frequency scale w 



ren considering 



the profile of I{uj) is UJ.^, near which this function attains 
its maximum. In turn this implies that 

\Ai{[[jr)\ « |A,(0)| - Huj,) - max{/(a;)}, (F18) 

and that in the interval (0, t^*), ^ (uj) > 0. On the other 
hand wc can use the mean value theorem to write 



7rA/(w) 





'dl 


/ dc^' 







duj 




f*OC 


Jo J 


3w/2 



-{n[u}', uj]) + 



1 



UJ + UJ 



-I{uj') 



UJ — Ul' UJ + U) 



-I{uj')duj' 



(F19) 



where ri[tj',w] e {uj,uj'). Equations ( |F17| -(|F19 1 imply 

(F20) 



|A/ {uj) I I{uj) ^ \Ai{Cjr) 



UJ UJ UjR 

for ujR < UJ < UJ■^,, and 

Ai (uj) « Ai (wh) , I{uj)<I (CbR) < I A/ {Cjr) \ , 

(F21) 

for UJ < lur. 

Furnished with relations ( |F17| ), ( |F20| ) and ( |F21| , we 
turn now to the analysis of defined in Equation 

(F15l. First we consider the case \uj — ujrI }^ ujr- For 
UJ > lor we use [w — (I'/jI ^ uj and ( |F20 ), while for uj < ujr 
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we use (F21 ). In both instances it is simple to establish 
|i?(a;)| < |A/ (ujh) |/w_r. To analyze the remaining case 
\uj ~ LUjil <^ Cofj we consider the Taylor expansions of 
the functions A/ (uj) and / (w) around the frequency Cor. 
Equation (F17l implies that substituting into Eq. (F15) 



the linear parts of these expansions results in a relative 
error for E(uj) of order \uj — Cjji\/LbR <ti 1. Thus we may 
write 



U! 



.di 



AliCoR) 



U) — ujR + i 



LOR + U) 



H^R.) 



(F22) 



which using again Eqs. (F17l and (F20l leads to the de 



sired result for E{uj). This completes the derivation of 
Eq. (82 1 for the frequency range in which it is used in 



Subsection VB It is worth noting that a straightforward 
derivation using Eqs. (|3T]), ([32]), ( |F11| ), and ( |F15[ ) leads 
to 



C(9)-C*(<?) 



E{uj) 



(F23) 



with 



^ , [(^g + ^r) i^^q - Wfl) + iURl{LOR)] 

2^LjRUJq 

(F24) 

which can prove useful when analyzing the manipulation 
of a specific resonance by coupling to a pseudospin^^. 

Finally, to warrant our derivation of Eqs. ([s]) and ([6| 
in Subsection |VB] we also need to establish that the func- 
tion I{uj) remains well behaved at w = ur as d/L — > 0. 
This can be accomplished by using the exact expression 
( 32 1 , Eq. ( 26 ) , and the behavior of the overlaps for small 



but finite d/L. To understand the latter it is essential 
to consider the qualitative behavior of corrections to the 
scattering modes This can be done"^ usmg: 

(i) that each of the refiection amplitudes admits an ex- 
pansion in powers of kd and those corresponding to a 
junction relate to the finite transmission coefficient via 
energy conservation and (ii) the exact reduction to a 
single junction performed in Sec. Ill We note that for 
our present purposes the precise form of the coefficients 
of the expansion (i), which will be studied elsewher^SS 



is not important. Propositions (i), {ii), and Eq. (26 1 al 



low us to write the following Lorentzian approximation 
L=l: 



S{L^)^ARid) 



TRid)/2 



uj~CjR{d)Y+Tl{d)/4: 
uj - ujRid) 



-LO~CoR{d)f+Tl{d)/4 



(F25) 



where the corrections are higher order in d for all fre- 
quencies in a small neighborhood of lor, and to establish 



lim AR{d) 



1 

2uj't 



lim BR{d) = 0. 



(F26) 



Inserting Eq. (F25l into Eq. (32 1 we obtain 



TR{d) 



LURARid) 2 



^BR{d) 



•^Rid)] 



Slid) 
Alid) 



[uj - u}R{d)] 



, (F27) 



which together with Eq. (F26l leads to I{uj) w TR^d) 



for frequencies close to lor implying that (as expected on 
physical grounds) the resonator mode's environment is 
structureless at the characteristic resonant frequency ujr 
(cf. Subsec.[!lB|. 

An analogous procedure allows to analyze the qualita- 
tive behavior at the other resonances lu = ujn =/= lur where 
we find that S^u) presents Fano profiles that result in cor- 
responding features in the function /(w) that have negli- 
gible relative spectral weight as c?/L —> 0. Thus, though 
the latter do not affect the behavior at other frequen- 
cies, they invalidate Eqs. (F20l, (F21l at these special 



frequencies precluding in their neighborhoods the use of 
Eq. (82l. Toconclude, we point out that from Eqs. (F25I, 
(84 1, and (80 1 one can obtain an approximation for the 



function SJlj) [in terms of TR{d)] adequate for all low 
frequencies other than uj„ ^ lor. The latter approxima- 
tion together with the asymptotic behavior discussed at 
the end of Subsec. VB and the exact relation (32 1 allow 



for an independent heuristic justification of Eq. (F17l 
— i.e. without resort to the approximation (82 1 or the 



properties of the spectral densities expected on physical 
grounds. 
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More precisely, at the origin taken as the center of mass 
of the beam's cross section at the junction. We note that 
the anomalous frequency scaling of the ppito) [cf. Eq. (58l] 
when viewed as the DOS of a ID fictitious environment 
expresses just the non-adiabatic nature of the junction. 
This yields kn,c/t = ■?i"(n + 5/2)/ L and k^^v/h ~ 7r('^ + 5 — 
l/2)/L where the relative error is of order e"*^" "/'*^. 
More precisely, they should yield the correct order of 
magnitude provided the corresponding resonance lies be- 
low the lowest infrared cutoff associated with the phonon 
waveguide afforded by the undercut — here we assume 
d(7 ~ i > t. 

The criterion underpinning the derivation is that the 
support-induced frequency shift (see below) be smaller 
than the "free spectral range" associated with the bare 
resonance. 

If there are support dimensions ds < d our derivations for 
the lowest order in d/ L remain valid provided the limit d —> 
yields a support satisfying assumption (i). Implicit in the 
latter is the idealization that for a bridge the two supports 
are effectively disconnected for phonons with wavelengths 
^ 27r/L. 

Note that the new set of modes {u'ii{f),u'g{'r)} is chosen 



(") 



= for im ~ 3 with m ^ 1. 



to be also orthonormal in the metric defined by Eq. ( |19[ ) 
— which accounts for the discrete mode's elastic energy 
density being proportional to ps rather than pt.. 
More precisely, its analog with ujji —^ LjR[ioii — A/(0)] 
[cf. Eqs. E^, EEl and Eq. ((83|]. This correction corre- 



sponds to the environment-induced renormalization of the 
resonator's frequency that does not cancel out given the 
lack of counterterm in Hamiltonian ( |20| l . 
We note that while u^(f)e~''^^' is only a solution of the 
elastic wave equation inside the resonator, the spatial aver- 
age of the elastic energy density should be taken over the 
whole structure when considering Hamiltonian H. How- 
ever, the boundary term that arises due to the disconti- 
nuity at the junction vanishes due to the choice of ho- 
mogeneous boundary conditions. The reader is referred to 
Ref. iSDI where a related treatment is given of the electro- 
magnetic analog of our problem. 

Note that we consider the limit d/L — + with k{q)L held 
constant. 

For the last two properties, in the case of the TRE solu- 
tions, we have also used that \ f\G{f — f' , u) is well behaved 
as ^ 0. 

These may be chosen to have well defined reflection sym- 
metries as for the free elastic half space. 
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This is just the analog for the decoupled support of Eq, 
and is obtained from it via the ansatz u ' , b — > 

The limits S ^ of the angles are taken keeping the shape 
of the cross section invariant. 

The 2D surface states (SAW) have a penetration depth 
that scales as 1/qt so that their contribution to the DOS 
has the same scaling with uj as the bulk contribution. 
We note that as the cross section S lies at 2 = one has 



Though the universal properties discussed in App. [bJ im- 
ply that for the end corrections the total force [zeroth or- 
der moment M'"'] and total torque [antisymmetric part of 
M'^-*] scale at most as the inertia, the symmetric part of 
M(i' and the M*^' can be zeroth order in k. 
There is naturally a purely real singular part that corre- 
sponds to the quasistatic near field and over an oscillation 
period does not perform net work. 

Note that the (f>n,i3{z) are homogeneous solutions of the 
linear systems defined by the RHS of Eqs. ( |34[ ), ( |36[ ) in the 
limit d/L 0. 

This disparity is just a trivial consequence of the disper- 
sion relations (which imply LJo,t/i^o,v/h ~ L/d) and the 
strong dependence of the corresponding Q-values with the 
aspect ratio — for similar dimensions the Q-values of the 
fundamental resonances of these branches are comparable. 
For asymmetric cross sections the (3 index is conventional. 
Any surface or edge states resulting from the free bound- 
aries have a penetration depth of order 1/qt- 
For a SWNT if we use the shell model the appropriate 
ansatz reads instead S — > 2-kR, Iz — > 2tv R{r'^)s — natu- 
rally, the results for the functions Pjsiuj) are independent 
of the "microscopic model" used [cf. Eqs. (47 1, (50l,(56l]. 



In instances where there are several disparate characteris- 
tic wavelengths — as is the case for the radial breathing 
mode of a microtoroid that hybridizes appreciably with 
the fiexural modes of the corresponding membran^^ — 
should be taken as the largest of the corresponding 
wavevectors. 

An example of particular relevance in the context 
of optomechanics is the radial breathing mode of a 
microtoroid^^'^. 

In the case of microtoroids with membrane thicknes^^^'^ h 
much smaller than the smallest pedestal radius one should 
instead define the resonator volume as the undercut region, 
so that S corresponds to the vertical rim at the smallest 
radius of the undercut. Thus, this scenario falls into the 



singular case (i) as defined below Eq. ( 10 1 and u'^ should 



instead be specified by clamped boundary conditions at S. 
It is straightforward to extend the derivation to resonances 
with other reflection symmetries R^, Ry — as for the sym- 
metric beam geometry each resonance only couples to the 
TRE branch with the same symmetries. 
To understand the normalization and the support's one- 
dimensional DOS (needed below to integrate over q) it is 
instructive to consider the infinite length limit of the elastic 
object generated by refiecting the pedestal with respect to 
S and identifying the resulting endpoints — with an ade- 
quate smoothening of the possible discontinuities in S' [z). 
Note that is normalized using the non-trivial metric 
defined by Eq. ([l9|. 



